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1.  INTRODUCTION  AND  EXECUTIVE  SUMMARY 


For  Ihe  period  from  15  May  1986  -  15  May  1988,  the  project  described 
in  this  report  sought  to  develop  a  more  complete  theoretical  frameuuork 
for  the  Plasma  Electron  Microuuave  Source  (PEMS)  concept  (1).  The  PEMS 
concept  is  a  novel  approach  to  generating  intense  pulses  of  microuuave 
power  by  transforming  energy  stored  in  a  mirror-confined  hot-electron 
plasma  into  large  amplitude  whistler  waves,  propagating  parallel  to  the 
static  magnetic  field  that  confines  the  hot-electron  plasma.  PEMS 
devices  can  function  either  as  amplifiers  of  externally  launched  whistler 
waves  or  as  oscillators,  depending  on  the  beta  and  temperature 
anisotropy  of  the  hot-electron  plasma,  where  beta  is  the  conventional 
ratio  of  plasma  kinetic  pressure  to  the  magnetostatic  pressure  of  the 
static  magnetic  field. 

The  main  ob jectives  of  the  theoretical  research  undertaken  in  this 
project  were  to  estimate  limiting  values  of  the  plasma  energy  that  could 
be  stored  efficiently  in  mirror-confined  hot-electron  plasmas  and  the 
fraction  of  that  stored  energy  that  could  be  transformed  into  RF  field 
energy.  The  method  of  hot-electron  plasma  generation  investigated 
here  is  electron  cyclotron  heating  (ECH)  including  advances  such  as  Upper 
Off-Resonant  Heating(UORH)  (2)  in  which  energetic  electrons  are  heated 
preferentially.  Alternative  methods  for  creating  the  hot-electron  plasma 
using  pairs  of  oppositely-directed  co-linear  steady-state  electron  beams 
are  also  under  active  investigation  at  AMPC.  Here  our  concern  is 
primarily  with  (1)  the  limiting  values  of  hot-electron  energy  density  set  by 
instabilities;  and  (2)  the  aspects  of  whistler-wave  growth  in  anisotropic, 
hot-electron  plasmas  affecting  the  gam  and  output  power  of  PEMS 
amplifiers  and  oscillators. 


The  report  consists  primariLy  of  tujo  technical  sections  describing 
results  achieved  in  each  of  these  tuuo  aspects  of  the  problem:  Section  2, 
"Formation  of  Stable,  High-Beta,  Relativistic-Electron  Plasma  Using 
Electron  Cyclotron  Heating"  (3),  has  been  published  recently  in  Nuclear 
Fusion.  Section  3,  "Rmplification  of  UJhistler  UJaves  Propagating  Through 
Inhomogeneous,  Anisotropic,  Mirror-Confined  Hot-Electron  Plasmas"  (41, 
has  been  published  in  Physics  of  Fluids.  The  most  recent  results 
achieved  in  this  project  uuhich  have  not  yet  been  submitted  for 
publication  are  reported  briefly  in  Section  4.  The  remainder  of  this 
section  provides  a  recapitulation  of  the  main  results  using  the  standard 
format  for  final  reports. 

A  one-dimensional,  steady  state,  relativistic  Fokker-Planck  model  of 
electron  cyclotron  heating  (ECH)  is  used  to  analyze  heating  kinetics 
underlying  the  formation  of  tuio-component,  hot-electron  plasmas 
characteristic  of  ECH  in  magnetic  mirror  configurations.  The  model  Is  first 
applied  to  uuell  diagnosed  plasmas  studied  experimentally  in  the  past  and 
then  used  to  simulate  the  efficient  generation  of  relativistic-electron 
plasmas  by  Upper  Off-Resonant  Heating  (UORHj  as  has  been  demonstrated 
empirically  (2).  The  characteristics  of  unstable  whistler  modes  and 
cyclotron  maser  modes  are  then  determined  numerically  for 
two-component  hot-electron  plasmas  sustained  by  UORH.  Cyclotron 
maser  modes  are  found  to  be  strongly  supressed  by  cold 
background-electron  species,  while  whistler  growth  rates  are  reduced  by 
relativistic  effects  to  levels  that  may  render  them  unobservable, 
provided  the  hot-electron  pressure  anisotropy  is  below  an 
energy-dependent  threshold.  A  fully  relativistic  local  dispersion  relation 
for  whistler  waves  is  solved  at  closely-spaced  points  along  the  magnetic 
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field  Lines  of  a  2:1  magnetic  mirror  in  luhich  anisotropic, 
spatialLy-inhomogeneous,  hot-electron  plasmas  are  confined.  The 
Limiting  plasma  parameters  for  convective  (spatial)  growth  are 
determined  numerically  and  used  to  identify  plasma  conditions  Leading  to 
maximum  amplification  of  microwave  signals  introduced  into  the  plasma 
in  the  form  of  whistler  waves.  The  maximum  gain  has  been  evaluated 
numerically  for  a  range  of  values  of  the  hot-electron  plasma  parameters 
within  which  all  known  stability  criteria  are  satisfied.  Very  high  gains  (up 
to  40  dB)  are  indicated  over  the  range  of  plasma  conditions  examined. 


1.1  RESEARCH  OBJECTIVES  OF  THE  PROJECT; 

(1)  Obtain  self-consistent  relativistic-electron  distribution  functions  for 
optimized  electron  cyclotron  heating  using  analytical  and  numerical 
solution  of  suitable  Fokker-Plenck  equations. 

(2)  Test  the  distribution  functions  obtained  in  (1)  for  violations  of 
equilibrium  and/or  stability  criteria  to  identify  plasma  processes 
limiting  the  maximum  achievable  stored  energy  density. 

(3)  Determine  the  effect  of  spatial  inhomogeneities  on  the  propagation 
and  growth  of  whistlers  in  the  interior  of  the  plasma. 

(4)  Evaluate  the  quasi-linear  behavior  of  unstable  whistler  waves  in  the 
interior  of  plasma. 

(5)  Model  the  wave  processes  governing  the  reflection,  transmission 
and  conversion  of  whistler  waves  at  the  plasma  surfaces. 
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1.2  ACCOMPLISHMENTS  AND  PROGRESS  TOUiARD  THE  OBJECTIVES 


(1)  Both  one-  and  tuuo-dimensionoL  Fokker-PLonck  models  of  relativistic 
ECH  have  been  solved  for  a  ojide  range  of  plasma  and  magnetic  field 
configurational  parameters,  including  toroidal  magnetic  geometries 
as  ujell  as  simple  operbended  magnetic  mirrors.  Fundamental  ECH 
and  Upper  Off  Resonant  Heating  results  from  experiments  in  ELMO 
and  SM-1  are  simulated  in  all  major  respects  by  the  Fokker-Planck 
model. 

(2)  Equilibrium  distributions  obtained  from  the  Fokker-Planck  models 
have  been  analyzed  in  detail  for  stability  against  electromagnetic 
modes  driven  by  pressure/temperature  anisotropy,  specifically  the 
so-called  "cyclotron  maser  mode"  and  unstable  uuhistler  uuaves.  In 
both  cases  only  uuaves  propagating  parallel  to  the  static  magnetic 
field  have  been  analyzed,  even  though  the  cyclotron  maser  mode 
achieves  its  maximum  grouuth  rate  for  normal  propagation.  It  is 
found  thot  these  modes  ore  completely  suppressed  by  the 
cold-electron  component  that  is  invariably  present  in  ECH  plasmas. 
The  uuhistler  uueve  grouuth  rate  is  governed  by  the  degree  of 
pressure  anisotropy,  the  hot-electron  beta,  and  the  degree  to  uuhich 
the  hot  electrons  are  affected  by  the  relativistic  increase  in  mass. 
Highly  relativistic  hot-electron  plasmas  are  found  to  exhibit  very 
small  uuhistler  grouuth  rates.  The  condition  for  existence  of 
equilibria  (uuithin  the  anisotropic  guiding  center  fluid  theory)  has  been 
evaluated  numerically  and  shouun  to  limit  the  degree  of  pressure 
anisotropy  more  stringently  thon  the  condition  for  absolutely 
unstable  uuhistler  uuaves  for  values  of  hot-electron  beta 
(perpendicular)  greater  than  15  -  20%. 

(3)  A  fully  relativistic  local  dispersion  relation  for  uuhistler  uuaves 
propagating  parallel  to  the  magnetic  field  has  been  solved  at 
closely-spaced  points  along  the  field  lines  of  a  2:1  magnetic  mirror  In 
uuhich  anisotropic,  spatiolly  inhomogeneous,  hot-electron  plasma  is 
confined.  The  limiting  plasma  parameters  for  convective  (spatial) 
growth  are  determined  numerically  and  used  to  identify  plasma 
conditions  leading  to  maximum  overall  amplification  of  microuuave 
signals  introduced  at  one  axial  edge  of  the  plasma  in  the  form  of 
uuhistler  uuaves.  The  maximum  gain  has  been  evaluated  numerically 
for  a  range  of  values  of  the  hot-electron  plasma  beta  and 
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temperature  anisotropy,  luithin  uihich  the  major  stability  criteria  are 
satisfied.  Overall  (global]  gains  of  40  dB  or  more  are  indicated  over 
the  entire  range  of  beta  values  investigated.  Banduiidths  in  excess 
of  15  -  20%  are  obtained  at  louj  beta;  and  much  larger  banduuidths 
are  found  in  the  more  relativistic  cases  luith  higher  betas. 

(4)  In  anticipation  of  rather  sophisticated  numerical  studies  of  the 
nonlinear  evolution  of  spatially  grouping  uuhistler  ujaves,  almost 
certainly  requiring  numerical  simulation  techniques,  lue  have  carried 
out  preliminary  estimates  of  the  maximum  energy  that  can  be 
transformed  into  luave  energy  using  thermodynamic  bounds 
developed  by  T.K.  Fouuler  (5],  Foujler  defined  a  type  of  free  energy 
that  is  appropriate  for  plasmas  capable  of  sustaining  microturt)ulent 
phenomena,  such  as  the  groujth  of  high-frequency  uuaves.  The 
magnitude  of  the  free  energy,  together  uuith  the  dynamical 
constraints  associated  uiiith  the  Vlasov  kinetic  equation  leads  to  a 
bound  on  the  maximum  amplitude  of  the  unstable  uuaves.  In  the 
present  case,  Fouuler's  estimate  for  the  limiting  amplitude  of  the  RF 
electric  field  energy  is  directly  proportional  to  the  hot-electron 
plasma  energy  and  the  degree  of  temperature  anisotropy  relative  to 
a  marginally  stable  case.  The  estimated  fraction  of  the  hot-electron 
stored  energy  that  can  be  transformed  Into  uuave  energy  for  the 
cases  examined  to  date  can  range  up  to  50%. 

(5)  UUe  have  identified  a  promising  spatially  varying  plasma 
configuration  that  offers  the  possibility  of  minimal  Internal 
reflection  of  uuhistler  uuaves  launched  at  one  end  of  the  plasma. 
The  hot-electron  plasma  is  confined  in  a  2:1  magnetic  mirror  uuith  an 
arbitrarily  long  flat  field  central  region.  The  effective  axial 
boundanes  of  the  hot-electron  plasma  are  taken  to  be  the 
fundamental  ECH  resonant  surfaces;  uuithin  this  central  region  cold 
plasma  is  generated  uniformly  through  ionization  of  incident  neutral 
gas.  The  cold  plasma  flouus  freely  along  the  lines  of  force  and 
expands  in  the  region  beyond  the  magnetic  mirrors,  uuhere  the 
magnetic  lines  of  force  diverge.  Since  the  propagation  velocity  of 
uuhistler  uuaves  approaches  the  speed  of  light  as  the  plasma  density 
decreases  by  virtue  of  the  diverging  magnetic  field  Lines,  uue 
anticipate  efficient  coupling  to  outgoing  vacuum  electromagnetic 
uuaves. 
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SECTION  2 


FORMATION  OF 

STABLE,  HIGH-BETA,  RELATIVISTIC-ELECTRON  PLASMAS 
USING  ELECTRON  CYCLOTRON  HEATING 

G.E.  GUEST,  R.L.  MILLER 

Applied  Microwave  Plasma  Concepts.  Inc., 

Carlsbad,  California, 

United  States  of  America 


ABSTRACT.  A  one-dimensional,  steady-stale,  relativistic  Fokker- Planck  model  of  electron  cyclotron  heating 
(ECH)  is  used  to  analyse  the  heating  kinetics  underlying  the  formation  of  the  two-component  hot-cleciron  plasmas 
characteristic  of  ECH  in  magnetic  mirror  conFigurations.  The  model  is  Ftrsl  applied  to  the  well  diaennscd  plasmas 
obtained  in  SM-I  and  is  then  used  to  simulate  the  effective  generation  of  relativistic  electrons  by  upper  olf  rcsonant 
healing  (UORH),  as  demonstrated  empirically  in  ELMO  The  characteristics  of  unstable  whistler  modes  and  cyelotion 
maser  modes  are  then  determined  for  two-component  hol-eleciron  plasmas  sustained  by  UORH  Cyclotron  maser 
modes  are  shown  to  be  strongly  suppressed  by  the  colder  background  electron  species,  while  the  growth  rates  of 
whistler  modes  are  reduced  by  relativistic  effects  to  levels  that  may  render  them  unobservable,  provided  the 
hot-electron  oressure  anisotropy  is  below  an  energy  dependent  threshold 


1.  INTRODUCTION 

Early  experiments  on  electron  cyclotron  heating 
(ECH)  of  plasmas  confined  in  magnetic  mirrors  (IJ 
using  CW  microwave  power  at  a  single  frequency 
created  an  unusual  plasma  containing  two  distinct 
populations  of  electrons.  The  majority  of  the  electrons 
remained  at  a  temperature  that  was  typically  below 
100  eV,  while  a  small  fraction  was  heated  to  average 
energies  above  100  keV.  Rather  broad  regimes  of 
stable  operation  were  found  to  be  limited  by  the  onset 
of  several  dominant  modes  of  instability,  such  as  the 
mirror  mode,  the  flute  mode  and  unstable  whistler 
waves.  Later  experiments  in  the  ELMO  facility  [21 
showed  that  the  energy  stored  in  the  hot-electron 
component  could  be  greatly  increased  by  using  upper 
off-resonant  heating  (UORH),  in  which  microwave 
power  at  a  frequency  above  the  cold-electron  gyro- 
frequency  was  added  to  the  fundamental  ECH  power. 
In  addition  to  increasing  the  relative  density  of  the  hot 
electrons,  UORH  also  enhanced  the  stability  of  the 
plasma,  even  for  values  of  beta  in  excess  of  50^  (beta 
is  the  conventional  ratio  of  plasma  kinetic  pressure  to 
magnetostatic  pressure). 

The  stable,  high-beta,  hot-electron  annulus  revealed 
by  the  ELMO  experiments  led  to  the  successful  use  of 
these  plasmas  (‘ELMO  rings')  to  stabilize  a  scalar 
pressure  core  plasma  confined  in  the  ELMO  Bumpy 


Torus  device  (see.  lor  example,  Rel  |.^|);  this  stimu¬ 
lated  widespread  interest  in  the  fundamental  properties 
of  these  remarkable  plasmas  (sec.  lor  example. 

Ref.  [4)).  Recently  proposed  applications  (or  these  hot- 
electron  plasmas  include  stabilization  of  localized 
modes  in  heliacs  |.S|  and  tokamaks  jb).  as  well  •" 
energy  storage  for  subsequent  conversion  to  high- 
power  microwave  pulses. 

It  has  required  the  efforts  of  many  workers  to 
evolve  reasonably  satisfactory  theoretical  models  of  the 
steady  state  equilibrium  produced  in  ECH  cxpcrirtictus. 
The  basic  characteristics  of  the  interaction  of  individual 
electrons  with  microwave  electric  fields  were  derived 
very  early  (for  a  discussion  of  early  theoretical  work, 
see  Ref.  |7))  and  useful  formulas  were  obtained  for  the 
resulting  heating  rate  under  a  wide  range  of  conditions 
(for  a  description  of  later  work,  see  Ref.  |8j>.  Sub- 
.sequcnl  development  of  geometrical  optics  codes  (sec. 
for  example.  Ref.  |9|)  made  it  possible  to  describe  the 
propagation  and  absorption  of  microwaves  in  specified 
plasma  media.  Progress  was  accelerated  by  the  formula¬ 
tion  of  a  relativistic  Fokkcr-Planck  theory  o(  ECH  by 
Bernstein  and  Baxter  jlO)  Hamasaki  and  co-vwirkcrs 
111)  took  momenta  of  the  Fokker  Planck  equation, 
using  trial  distribution  function^,  in  an  attempt  to 
predict  the  power  required  for  fointation  of  the 
ELMO  rings.  Other  workcis  |I2|  developed  numerical 
tcchtiiqucs  for  a  direct  stviution  ('f  the  Fokker  I’l.inck 
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equation;  these  are  now  implemented  in  large  computer 
codes  [13).  Ultimately,  these  codes  are  expected  to 
provide  a  detailed,  self-consistent  picture  of  high-beta 
hot-electron  plasmas  formed  by  FX'H, 

Nonetheless,  there  is  an  immediate  need  for  a 
reduced  model  that  can  be  used  to  search  over  the 
wide  range  of  parameters  characterizing  an  ECH 
experiment.  It  is  not  practical  to  carry  out  such  a 
search  with  two-dimensional  Fokker- Planck  codes; 
this  would  require  broad,  time  consuming  parameter 
searches  to  overcome  the  general  lack  of  experimental 
measurements  of  such  essential  parameters  as  the 
strength  of  each  RF  electric  field  component  present  in 
the  plasma.  We  have  addressed  this  need  by  evolving  a 
Fokker-Planck  model  for  the  steady  slate  distribution 
of  electrons  in  energy;  this  model  includes  in  a  heuristic 
way  the  phenomena  that  depend  essentially  on  the 
electron  orientation  in  a  gyrotropic.  two-dimensional 
velocity  space.  This  model  has  been  described  in 
considerable  detail  in  Ref.  (14). 

In  order  to  model  a  complex  physical  situation  in  a 
relatively  simple  way,  the  power  balance  issues  affecting 
the  formation  of  a  relativistic  electron  tail  have  been 
considered  primarily  and  not  the  particle  balance 
issues,  except  in  so  far  as  a  significant  associated 
energy  loss  occurs,  as  in  the  case  of  non-adiabatic 
scattering  of  high  energy  electrons  into  the  loss  cone. 
For  a  more  detailed  justification  of  this  approach,  sec 
the  Apnendix.  Thus,  we  have  chosen  to  c(>nserve 
particles  while  accounting  for  the  energy  loss  by 
converting  the  non-adiabatic  loss  term  into  a  dynamical 
friction  term,  as  described  in  Appendix  C  of  Ref.  (14). 
This  choice  has  the  desired  physical  effect  of  sharply 
reducing  the  hol-eicciron  distribution  function  for 
energies  above  the  threshold  for  violation  of  adiabatic 
invariance,  as  is  clearly  evident  in  Figs  la  and  lb.  The 
hot-electron  distribution  in  energy  is  shown  in  Fig.  la 
for  plasmas  typical  of  those  reported  in  Ref.  [15j. 
Figure  lb  shows  the  four  rates  accounted  for  in  the 
model;  the  quasi-linear  heating  rate,  the  Coulomb 
thermalization  rate,  the  synchrotron  radiation  cooling 
rate  and  the  energy  loss  rate  due  to  non-adiabatic 
scattering. 

Here,  we  first  apply  this  reduced  model  to  the  SM-1 
experiments  [l.S|  in  order  to  demonstrate  the  success  of 
this  model  in  describing  the  hot-eiectron  plasmas 
studied  in  SM-1.  The  results  of  this  demonstration  are 
presented  and  discussed  in  Section  2.  They  clarify  a 
number  of  features  of  ECH  experiments  that  had  not 
been  fully  understood  previously.  In  particular,  the 
role  of  microwave  electric  fields  parallel  to  the  static 
magnetic  field  (ordinary  mode  of  propagation)  is 
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SMuhriiirnn  radtanon,  and  rnrriiv  h\  non-ailitihaiti  xiatwrine 
of  hot  cicctrnns  into  ifw  lox  x  cone. 


shown  to  be  critical  to  an  cffectise  formation  of  the 
hot-electron  component.  Furthermore,  by  exploiting 
the  scnsilixily  of  the  hot  electimi  heating  rate  to  the 
parallel  microwave  electric  fields,  we  are  able  to  use 
the  reduced  Fokker- Planck  FX'H  model  to  demonstrate 
the  impact  of  non-adiabatic  scattciing  [Ibl.  which  has 
been  conjectured  lor  some  time  to  limit  the  maximum 
temperature  of  the  hot  electrons  in  most  ECH  experi¬ 
ments  to  date. 
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We  then  apply  the  present  model  to  the  UORH 
experiments  in  ELMO  (2|,  where  the  model  displays 
the  same  striking  increases  in  stored  energy  with 
increasing  UORH  power  that  were  observed  in  ELMO. 
Finally,  the  model  is  used  to  obtain  high-beta  hot- 
electron  equilibria  suitable  for  subsequent  stability 
studies.  Although  the  reduced  model  cannot  provide 
information  on  pressure  anisotropies  or  other  two- 
dimensional  properties  of  the  equilibria,  it  does  relate 
the  relative  temperatures  and  densities  of  the  two 
electron  species  and  provides  a  much  larger  degree  of 
internal  consistency  than  was  heretofore  available. 

In  Section  3,  we  analyse  the  stability  of  the  hot- 
electron  plasmas  to  transverse  waves  propagating  along 
the  static  magnetic  field,  namely  whistler  waves  (171 
and  the  cyclotron  maser  mtxies  (see,  for  example. 

Ref.  (18a]  and  the  excellent  set  of  references  cited 
therein;  for  more  recent  results,  see  Ref.  (18b|).  It  has 
been  known  for  many  years  (19[  that  relativistic  effects 
generally  reduce  the  growth  rates  of  unstable  whistler 
waves,  but  specific  results  depend  on  the  assumed 
equilibrium.  Here  we  u.se  the  Fokker-Planck  model  to 
obtain  equilibria  that  aie  consistent  with  ECH  kinetics 
and  demonstrate  the  relativistic  suppression  of  whistler 
growth  rates.  The  cyclotron  maser  modes  arc  found  to 
be  completely  stabili/ed  by  the  cold-electron 
component. 

Finally,  we  note  in  passing  that  our  understanding 
of  flute  modes  in  high-beta  hot-electron  plasmas  has 
not  been  considered  here  and  i  "mains  incc'inplcte. 
although  fully  electromagnetic  kinetic  models  of  line 
tying  (20[  have  helped  clarify  finite-length  effects 
There  is  an  active  effort  n  predict  the  stability  of  these 
high-beta  plasmas  to  flute-bkc  modes  (21  (  which  wc 
hope  will  he  aided  by  more  detailed  descriptions  of 
realistic  equilibria.  Discussion  and  conclusions  arc 
found  in  Section  4. 

2.  FORMATION  OF 

HOT-ELECTRON  PLASMAS  USING  ECH 

In  this  section  we  describe  the  results  ol  applying 
the  reduced  Fokker-Planck  ECH  model  to  liCH 
experiments  in  the  mirror  devices  SM-1  (l.‘'l  and 
ELMO  (2)  before  using  the  model  to  obtain  high-beta 
equilibria  for  stability  analysis  The  Fokker-Planck 
model  is  described  in  great  detail  in  Ref.  |I4|  and 
references  cited  there 

Briefly,  wc  use  a  nnc-dimcnsion:il.  steady  stale, 
bounce-averaged  Fokker-Planck  rrtodcl  of  the  FCH 
process.  Since  our  primary  concern  is  with  the  llicrmali- 


TABLE  I.  SM-1  RFFERFNCF  CASK  (Ref  (l.'^li 


SlnrcJ  cncrcs 

W 

--  :n  .1 

R:ufttis  of  annulus 

f , 

(n  i  cm 

Radial  v^idth 

= 

^58  cm 

Arial  length 

1. : 

=  12  and  18  cm 

Annulu';  volume 

V..I 

i:  •  in'  cm  ' 

Mirror  ratio  on  axis 

M,,  - 

-  ■>  *> 

Coil  separation 

I,  = 

71  cm 

Mirror  ratio  at  r,i 

M  - 

Mirror  ratio  at  resonance 

M,  , 

—  1 

ECH  power 

P  == 

1  kW 

Micr<»wavc  frequenev 

f  = 

y  .S  (-.11/ 

Bat  kpr<nind  plasma 

•Jcnsiiv 

ft 

[  1  1  '  111 ' '  cm 

Icmpcralure 

1  = 

lo-.sn  cV 

electron  compimeni 

vtcnsilv 

n,  « 

-  10  '  (.Ill  ’ 

temperature 

b.  = 

.v.so.ann  keV 

ration  of  a  relatively  small  concentiaiion  id  energetic 
clcclri'iis  we  adopt  a  lincarired.  lest  pailicic  Iraiiie 
work  I  he  relativistic  Choii  eleiiioiis  aic  then 
regarded  as  a  distinct  species  whose  spatial  density  n,, 
is  much  smaller  than  the  dcnsltic^  o(  the  other  plasma 
species  All  ol  these  'fieUr  specie  •  are  assuiiicd  tv>  be 
isotiopic.  thermal,  and  at  a  much  lowei  lempeiatme 
than  the  hot  electrons  This  model  determines  the 
resp(>nsc  of  a  small  relative-density,  hirh  enerpv 
electron  tail  to  the  incident  miciowave  power  bv 
including  the  dominant  dynamical  processes  quasi 
linear  healing,  ibei  maliration  by  (.'(nilomb  scattering  (>n 
cooler  cleclrfvns,  synchrotron  radiation,  and  ilirect  loss 
by  iK'n  adiabatic  scaiteiing  Specilically ,  the  model 
predicts  the  electron  distribution  in  energy  resulting 
from  RF-  electric  fields  with  spccitied  v.ilucs  of 
frequency,  amplitude  and  polan/ation.  mciileni  upon  a 
baikprouiitl  plasma  ol  spci  ilied  di  iisitv  and  tempeia 
line,  ci’iilined  in  a  presctd>ed  mil roi -like  m.ignelis 
field 

\k  e  lake  as  iwir  SM  I  eleience  v.ise  ih.-  2('  I 
equilibrium  discussed  at  lengili  m  Rel  !1.‘'|  and 
summari/ed  in  I  able  I  .As  menii''iud  e.iiliei.  iheie  aic 
no  direct  mcasuiements  dl  the  siiciigih  I'l  the  imcio 
w.ne  electric  fields,  and  we  iiuol  ^laii  (wer  .ill 
plausible  values  ol  both  component--  I  ,ind  I  Vk  c  ilo 
ihi-  bv  spci  living  the  r.iiio  I-,  I  ol  the  ()  moile  and 
X-niode  lield  stiengths  and  seats Imig  o-n  the  values  o| 
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FIG.  2n.  Perpendicular  component  of  the  nn<  ronacc  electric  field, 
h  ,  ,  averaee  energy  of  the  hot  electron'..  T,,.  and  total  RF  power 
density.  Pf,.  versus  background  plasma  temperature.  1].  for  an 
assitmed  ratto  of  0-mode  to  X-inodc  field  strengths.  E,  F.  ^  =  2. 
tn  the  SM-I  configuration  summarized  in  Table  I  The  relative 
hot-electron  densitx  ts  n,,  n  =  0.!. 


FIG  2h  Variation  of  E  .  7"^  and  RF  pout  density  Ppf  uiih 
F,  F.^  for  a  background  plasma  temperature  T,  =  .10  eV  and 
n^  n  =  0.1 


required  to  sustain  the  observed  relati\e  hot- 
electron  density,  n^/n  =  0.1;  this  yields  the  observed 
hot-electron  average  energs  and  the  RF-  pc'vscr 
density  Pr,  required  to  sustain  the  equilibriiiin.  Since 
the  microwave  power  is  launched  piiniarily  in  ordinarv 
modes  and  is  multipis  reflected  from  the  casity  walls, 
w'c  assume  an  essentially  isotropic  tlistribuiK'n  of  wave 


vectors  k  for  the  ECH  wa\es.  An  illustialisc  case  is 
shown  in  Fig,  2a.  which  displays  F_ .  T,,  and  fVi 
versus  the  background  electron  tcilipcraturc  I  ,  for 
E|/E.  =  2.  The  chtiractcristic  scale  length  of  the  static 
magnetic  field  has  been  set  at  1.,  -  .^0  cni  for  this 
case,  a  choice  that  we  will  discuss  at  length  in  the 
following.  Nt'te  that  the  RF  power  density  required  to 
sustain  this  equilibrium  is  around  20  mW  em  '. 
accounting  for  roughly  250  W  of  RF  power.  This  docs 
not  include  the  fundamental  ECH  power  necessary  to 
sustain  the  background  plasma.  If  the  energy  e.xpended 
in  creating  each  ion-electron  pair  is  roughly  50  eV .  we 
estimate  this  fundamental  ECH  power  requirement  at 
25  mW  em"’  for  an  ambient  hydrogen  gas  pressure  of 
5  X  10“*’  torr.  Since  the  cold-plasma  volume  is 
estimated  at  25  L  in  Ref.  115).  as  much  as  625  \V  of 
fundamental  ECH  power  may  be  required  simply  to 
sustain  the  background  plasma  and  only  250  W  may  be 
required  to  sustain  the  hot-electron  component.  The 
variation  of  E^.  T|,  and  with  the  ratio  E|/E,  is 
shown  in  Fig.  2b  for  a  background  electron  tempera¬ 
ture  =  30  eV.  Note  the  rapid  rise  in  RF  power 
density  for  E,/Ej  >  10.  in  which  Ppp  increases  five 
times  while  the  average  energy  of  the  hot  electrons 
remains  almost  constant  The  saturation  of  despite  a 
substantial  increase  in  RF  power  indicates  the  rapid 
increase  in  the  non-;nliabatic  scaiieiing  rtite.  II  these 
calculatitms  arc  repeated  fim  difletcnt  \alucs  of  the 
magnetic  scale  length  Lf,  we  obtain  the  striking  results 
shown  in  Fig.  ?a.  It  is  especially  instructive  to  display 
the  Li  dependence  of  the  saturated  value  of  the  hot- 
electron  average  energy  and  the  corresponding  values 
of  the  electron  gyroradius: 

p.„  =  1.7  X  lO'  Cl  cm  \  (t2  -r  fl'B 

where  (  =  Th.,.'m<,c-  and  B  is  the  magnetic  intensity  . 
Since  we  require  p  £  kLi  for  adiabatic  confinement  of 
electrons,  we  can  evaluate  k  by  plotting  versus  L.. 
The  value  of  k  obtained  from  Fig.  3b  is  0.06.  which 
agrees  reasonably  well  with  the  ctmventional  result  cited 
by  Hckan  |I6|.  From  the  values  obttiincd  for  I,,  ...  it 
appears  that  Li  =  24  cm  is  a  reasonable  value  of  the 
magnetic  scale  length  in  the  hi't-clcctron  region  of 
SM-I  From  a  partibolie  fit  to  the  vacuum  magnetic 
intensity  we  wcvuhl  estimate  L>  =  21  cm.  in  reasonable 
agreement.  Nt'ic  that  this  is  much  less  than  R  ,  the 
radius  of  curvtiturc  of  the  magnetic  field  lines  in  this 
region,  which  is  estimated  to  be  .irouiul  bi'  cm 
I'or  large  values  of  I  particularly  at  higher 
magnetic  fields,  synchrotron  i.idialimi  becomes  .i 
dominant  loss  mechanism  in  the  model  fhis  h'ss  is 
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FIG.  3a.  Tj  Vfrsus  P^f  for  different  values  of  the  maftnetic  scale 
length  Ly. 


FIG  3h  Saturated  \alue  of  T,,  and  rorrespnndins  value  ol  the 
average  hot-electron  g\roradiiii  p,„,  versus  Li  The  wild  line  is  for 
P,„  -  0  06  L, 


circumvenied  in  the  case  discussed  hete  b\  spcci(\ing 
values  of  L|  similar  to  the  ELMO,  (high  field)  experi¬ 
ments;  L|  =  10  cm  (ELMO). 

Another  important  feature  of  the  dependence  of  T,, 
on  E|/E^  is  the  rapid  initial  increse  in  7'^.  with  F'r, 
constant  or  even  decreasing.  This  behaviour  makes 
clear  the  substantial  benefit  of  increased  microwave 
integrity  of  the  vacuum  chamber.  If  E|/Ej^  >  10  can 
be  sustained  in  a  suitably  designed  cavity  where  wall 
reflections  do  not  rapidly  convert  ordinary  modes 
to  extraordinary  modes,  can  be  increased  by 
roughly  100  keV  or  more.  Since  the  value  of  E, 
required  to  achieve  this  level  of  heating  is  approxi¬ 
mately  10^  V  cm'',  while  the  power  density  dissi¬ 
pated  in  the  plasma  is  around  1  W  cm'\  the  effective 
quality  factor  must  exceed  300.  The  ELMO  device  [2] 
utilized  a  very  high-Q  cavity  together  with  UORH 
power  to  obtain  exceptionally  efficient  heating  of 
relativistic  electrons.  The  frequency  of  the  UORH 
power  is  chosen  to  be  above  the  cold  electron  gyro- 
frequency  throughout  the  entire  cavity.  This  power 
is  therefore  absorbed  primarily  by  the  relativistic 
electrons  and  only  negligibly  by  the  background 
electrons.  Even  for  the  hot  electrons,  the  absorption 
rate  is  small;  thus,  if  the  cavity  has  a  large  quality 
factor,  the  electric  field  strength  can  increase  to  large 
enough  values  to  render  the  heating  rate  (and  hence  the 
absorption  rate)  appreciable.  Since  heating  rates  are 
quadratic  in  electric  field  strength,  they  can  be  large 
enough  to  sustain  high  energy  electron  populations  at 
levels  below  KX)  V  cm”'. 

The  Eokker-Planck  simulation  of  early  ELMO 
experiments  on  UORH  |21  are  shciwn  in  Fig.  4.  Here 
the  fundamental  ECH  is  at  a  frequency  of  10.6  GHz 
while  the  UORH  source  has  a  frequency  of  ^5.1  GHz. 
These  results  compare  favourably  with  the  empirical 
results  displayed  in  Fig.  4  of  Ref.  [2).  The  background 
plasma  density  and  temperature  are  2  x  10"  cm  ' 
and  50  eV,  respectively,  and  we  have  arbitrarily  set 
E|/Ei  =  10  for  both  sources.  In  the  absence  of 
UORH  power,  the  relative  density  of  hot  electrons  is 
roughly  2.5%.  sustained  by  an  RE  power  density  of 
less  than  6  mW-cni  '.  This  relative  density  increases 
to  64%  and  the  hot-electron  beta  approaches  50%  as 
the  UORH  power  is  increased. 

Having  calibrated  our  model  on  .SM-I  and  ELMO, 
it  is  of  interest  to  see  how  much  hot-electron  energy 
can  be  stored  in  comparably  sized  devices  by  optimizing 
the  parameters.  By  using  UORH  at  higher  frequencies 
(60'9()  GHz)  in  a  device  that  is  similar  in  physical 
scale  to  ELMO,  we  can  obtain  Eokkct-I'lanck 
equilibria  with  average  hot-eicction  cncigics  above 
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FIG.  4.  Average  hot-elearon  energy,  T^.  relative  hot-electron  den- 
sin.  ni,/n.  and  plasma  beta.  0.  versus  power  density  from  the 
UORH  source.  PmitH-  parameters  desrrihing  the  background 
plasma  and  the  heating  conftgitration  are  given  in  the  text. 


nhc^.  The  relativistic  effects  on  plasma  stability  can  be 
clearly  examined  in  this  hot-electron  temperature 
range.  In  these  studies,  we  have  again  specified  the 
background  plasma  density  at  n  =  I  x  10”  cm  ’  and 
varied  the  background  electron  temperature  from 
250  to  3000  eV.  We  assume  that  optimal  background 
electron  temperatures  can  be  maintained  by  a  suitable 
level  of  fundamental  ECH  power,  depending  mainly  on 
the  confinement  properties  of  the  configuration.  For 
the  Fokker-Planck  model  of  the  hot-electron  population 
we  reduce  the  fundamental  ECH  electric  field  strength 
to  a  negligible  level:  10"’  times  the  UORH  field 
strength.  We  then  vary  T,.  over  its  range  to  determine 
the  values  of  UORH  power  density  Prf  required  to 
sustain  specified  relative  densities  of  hot  electrons.  The 
Fokker-Planck  model  evaluates  T^  for  specified  values 
of  the  various  system  parameters,  such  as  the  magnetic 
field  scale  length  L>.  here  set  at  10  cm,  as  in  ELMO. 
These  results  are  displayed  in  Fig.  5.  Note  that,  since 
these  curves  are  deduced  for  constant  hot-electron 
densities  Ah  =  nn/n.  the  resulting  value  of  beta  is 
simply  porportional  to  Th-  The  RF  power  density 
exhibits  a  broad  minimum  that  moves  from  500  eV  to 
800  eV  as  the  hot-electron  density  ratio  is  increased 
from  10%  to  40%.  The  associated  beta  values 
approach  50%  and  the  hot-electron  average  energy  is 
around  550  keV.  These  are  the  equilibria  whose 
stability  will  be  studied  in  the  following  sections. 


3.  SIABILIIY  ANALYSIS 
FOR  TRANSVERSE  WAVES 
PROPAGATING  PARALLEL  TO 
THE  STATIC  MAGNETIC  El  ELD 

As  discussed  in  Section  2,  the  first  task  addressed  in 
this  work  was  to  obtain  high-beta  equilibria  from  the 
reduced  one-dimensional  Fokker-Planck  model  of  the 
ECH  process.  Guided  by  the  empirical  results  of 
Dandl  et  al.  [2]  on  UORH  in  the  ELMO  experiments, 
we  found  a  corresponding  regime  in  the  Fokker-Planck 
analysis  that  yields  high-beta,  relativistic-electron 
equilibria  ideally  suited  to  a  stability  analysis,  namely 
beta  values  around  50%  and  relativistic-electron 
temperatures  above  500  keV.  Our  one-dimensional 
Fokker-Planck  model  cannot  provide  any  information 
on  the  temperature  anisotropy  that  is  characteristic  of 
hot-electron  plasmas  formed  by  ECH  and  which  plays 
a  dominant  role  in  the  stability  of  these  plasmas  |22]. 
Fully  two-dimensional  Fokker-Planck  ECH  codes 
are  now  available  and  will  eventually  provide  two- 
dimensional  equilibria.  These  numerical  equilibria  can 
be  expressed  analytically  as  a  superposition  of  the 
equilibrium  distribution  functions  described  by  Ciiicst 
and  Dory  123): 
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where  f  =  0,  I,  2.  and  u  =  p/m  is  the  momen¬ 
tum  per  unit  rest  mass.  Here,  we  examine  our  one¬ 
dimensional  equilibria  for  stability  against  transverse 
waves  propagating  parallel  to  the  static  magnetic  field, 
with  the  temperature  anisotropy  rj  =  Tj^  /T|  and  the 
loss-cone  index  f  set  at  arbitrary  values.  We  single  out 
the  whistler  mode  [17,  19,  24)  and  the  cyclotron  maser 
mode  [18,  24]  because  of  the  relative  stringency  of  the 
stability  criteria  for  electrostatic  and  whistler  modes,  as 
discussed  by  Guest  and  Sigmar  [22],  and  because  the 
cyclotron  maser  mode  was  omitted  from  consideration 
in  Ref.  [22]. 

In  what  follows,  we  describe  an  efficient  and 
very  rapid  numerical  technique  for  solving  the  fully 
relativistic  dispersion  relation  (or  right-hand,  circularly 
polarized,  transverse  waves  propagating  parallel  to  the 
equilibrium  magnetic  field.  The  derivation  of  the 
dispersion  relation  is  sketched  here  for  completeness; 
more  detailed  discussions  arc  readily  available  [25]. 

The  relativistic  Vlasov  equation  for  electrons  is 


X  u,(E,e  -t-  E.e"') 

where  <i>  is  the  polar  angle.  In  the  uniform  magnetic 
field  assumed  here,  the  unperturbed  electron  orbits  are 
given  by  z  =  Z(|  -I-  Upt/y  and  </>  =  (i>o  +  eBt/ym,  with 
U|  and  Uj^  constant.  Choosing  Zq  =  0  at  t  =  0,  we 
obtain 
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From  Faraday  's  and  Ampere's  laws. 

V  X  (V  X  E)  =  -no(d/dl)ir  +  fn(dE/dt)] 


(dfidl)  +  (u/yltdf/dr) 


-  (e/m)[E  (u  X  B)/>](3f/du)  =  0 


k-E  -  (k-E)k  -  (w-/c‘)E  =  iu'^J 


where  =  1  -(-  u’/c’  is  the  Lorentz  factor  and  the 
velocity  is  v  =  u/y.  The  perturbed  distribution 
function,  associated  with  a  wave  whose  electric  and 
magnetic  fields  are  E  and  B  =  k  x  E/w,  is  given  by 

f|{r,  u,  t)  =  ^  1  dt' 
m  J-« 

^  \  -I-  k- 

\  ys>)  /  du  yw  du 


For  purely  transverse  waves  propagating  parallel 
to  the  equilibrium  magnetic  field,  with  components 
E,  =  (E,  ±  iEj/v'2.  this  expression  becomes 


Since  k-E  =  0  for  transverse  waves,  E  and  j  must 
satisfy 

l(k'c-/a!‘)  -  1)E  =  (i/u.'<n)j 

For  the  right-hand,  circularly  polarized  component,  the 
resulting  dispersion  relation  is 
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For  the  class  of  equilibrium  distribution  functions 
discussed  by  Guest  and  Dory  (23),  the  numerator  of 
the  integrand  can  be  rewritten  conveniently  as 


Since  the  Lorentz  factor  y  appears  in  the  resonant 
denominator  and  since  we  are  concerned  with  strongly 
relativistic  plasmas  with  y  s  2.  it  is  necessary  to 
evaluate  the  integral  numerically.  We  now  consider 
integrals  of  the  general  form 


l(u,  k)  =  1  duj  dui  e 


L  2 


F(U|.uj  .oj.k) 
kui  -  70!  +  eB/m 


For  complex  values  of  w  and  k,  there  may  be  zero, 
one  or  two  solutions  to  the  relativistic  resonance 
condition  kui  -  701  +  eB/m.  which  we  denote 
by  U|,  (i  =  1,  2).  We  denote  each  solution  for 
U|(o),k,ui)  by  w  and  the  corresponding  values  of  7 
by  7?  =  V 1  +  w‘/c^  +  u^/c^  Now  consider  the 
following  integral: 


where  Z  is  the  conventional  plasma  dispersion  function 
[26).  From  this,  we  obtain: 
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With  the  singularities  removed  in  this  manner,  the  U| 
integral  can  be  evaluated  using  the  Gauss-Hermite 
quadrature  [27): 

I  du,  E(u,)  =  W'h  ,.  E(u,,,) 

J  -  »  ,  =  i 


Similarly,  the  u^  integral  can  be  evaluated  using  the 
Gau,ss-Laguerre  quadrature  [26j 


dui  e'“^''“'’l'(uf)  =  W,  ,'F(u;  ,i,t 


To  evaluate  the  dispersion  relation  for  given  values  of 
w  and  k,  we  then  choose  a  value  of  s  (s  =  12  typically) 
and  for  each  of  the  s  prescribed  values  of  u;  , 

(i  =  1 . s)  we  evaluate  the  corresponding  root  or 

roots  of  the  relativistic  resonance  condition,  u,  ,  =  w,. 
such  that 

+  cB/m  =  0 

The  Gauss-Hermite  formula  is  then  used  to  calculate 
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FIG.  6.  Density  dependence  of  cyclotron  maser  mode  growth 
rate  for  fixed  values  of  beta.  The  plasma  parameters  are  specified 
in  the  text. 


A  Standard  routine  is  used  to  evaluate  the  plasma 
dispersion  function,  and  the  Gauss-Laguerre  quadrature 
is  used  to  evaluate  the  uj  integral; 


like/loss-cone  aspect  of  the  mirror-confined  relativistic- 
electron  equilibria  [28).  Of  particular  interest  for  the 
present  studies  is  the  stabilizing  effect  of  a  cold- 
electron  component  on  the  cyclotron  maser  mode. 

Since  our  Fokker-Planck  ECH  equilibria  consist  of  a 
minority  population  of  relativistic  electrons  in  a 
majority  population  of  much  colder  background 
electrons,  cyclotron  maser  modes  are  anticipated  to  be 
strongly  suppressed.  This  is  illustrated  in  the  results 
shown  in  Fig.  6,  where  we  fix  the  beta  value  at  50% 
and  100%,  and  the  relative  density  of  hot  electrons 
at  30%.  The  hot-electron  distribution  has  a  loss 
cone  index  f  =  2.  and  a  temperature  anisotropy 
jj  =  T^/Ti  =  3.  The  peak  temporal  growth  rate  W| 
normalized  by  the  (non-relativistic)  electron  gyro- 
frequency  a)„  is  plotted  against  the  total  electron 
density,  specified  through  the  ratio  of  electron  plasma 
and  gyrofrequencies,  For  the  parameters  of 

interest,  the  peak  growth  rates  are  less  than  10“’  times 
the  gyrofrequency.  Figures  7a  and  7b  demonstrate  the 
fundamental  dependence  of  the  cyclotron  maser  modes 
on  the  loss-cone  index  f,  rather  than  the  temperature 
anisotropy.  For  these  cases,  the  hot-electron  beta  is 
fixed  at  30%,  the  total  density  at  =  Wf,/2,  and  the 
relative  density  of  hot  electrons  at  50%. 
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F/G.  7a.  Demonstration  of  the  independence  of  the  cyclotron 
maser  mode  growth  rate  with  respect  to  temperature  anisotropy. 


With  this  algorithm  for  evaluating  the  velocity-space 
integrals,  we  can  use  standard  root-finding  routines  to 
solve  the  dispersion  relation  for  fully  relativistic  hot- 
electron  plasmas.  In  this  way  we  can  evaluate  the 
threshold  conditions  and  growth  rates  for  unstable 
whistler  and  cyclotron  maser  modes  in  the  high-beta 
hot-electron  plasmas  obtained  from  our  Fokker-Planck 
model.  As  has  been  pointed  out  by  Chu  and  Hirshfield 
[18a],  whistler  waves  with  phase  velocities  below  the 
speed  of  light  are  destabilized  by  the  temperature 
anisotropy,  while  the  cyclotron  maser  mode, 
corresponding  to  waves  with  phase  velocities  greater 
than  the  speed  of  light,  are  destabilized  by  the  beam- 


FIG.  7h  Dependence  of  the  cyclotron  nuner  mode  growth  rale  (>n 
the  loss-cone  indi  i  f 
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FIG.  8.  StabiUzation  of  the  cyclotron  maser  mode  by  background 
cold  electrons.  The  relative  density  of  hot  electrons  is  A*  =  n^/n. 


Finally.  Fig.  8  shows  how  ihe  growth  is  suppressed 
by  decreasing  the  relative  density  or  hot  electrons.  Ah. 
In  this  particular  case,  stabilization  is  complete  for 
Ah  S  40%;  here,  0  =  30%  and  -  0.5.  Tsang 

(24]  considered  this  effect,  but  only  for  large  values  of 
Ah  where  he  did  not  see  complete  stabilization. 

By  comparison,  the  whistler  waves  can  have  growth 
rates  that  are  larger  by  two  orders  of  magnitude  if  the 
anisotropy  exceeds  a  critical  value,  which  is  about 
three  for  parameters  characteristic  of  hot-electron 
plasmas  formed  by  ECH.  Indeed,  whistler  growth  rates 
can  be  reduced  to  low  values  mainly  by  increasing  the 
average  energy  of  the  hot  electrons  and  decreasing 
the  temperature  anisotropy.  These  two  effects  are 
illustrated  in  Figs  9a  and  9b,  in  a  single-species  hot- 
electron  plasma  with  f  =  0,  rj  =  3  and  0  =  30%. 

The  average  perpendicular  energy  is  expressed  as 

=  Th/mc^  in  Fig.  9b,  =  1.  Finally,  we  evaluate 

the  whistler  growth  rates  in  a  two-component 
equilibrium  obtained  from  the  Fokker-Planck  ECH 
model:  Th  =  527  keV,  An  =  30%,  T,  =  800  eV, 

=  0-5  and  0  =  31%.  The  temperature 
anisotropy  has  been  varied  (r;  =  3.  6.  10.  20)  to  show 
the  rapid  decrease  in  growth  rate  for  decreasing 
anisotropy,  as  displayed  in  Fig.  lOa.  For  the  weakly 
unstable  whistler  waves  the  growth  is  convective.  It  is 
thus  of  interest  to  evaluate  the  spatial  amplification 
rates  k,.  as  illustrated  in  Fig.  lOb. 


4.  DISCUSSION  AND  CONCLUSIONS 

We  have  demonstrated  several  points  in  this  paper. 
First  of  all.  using  the  classical  Fokker-Planck  miKlel  of 
ECH,  it  is  possible  to  reproduce  all  ol  the  features  of 
the  SM-1  experimental  results  that  can  be  compared 
with  the  reduced  one-dimensional  model.  In  particular, 
there  is  no  need  to  invoke  collective  processes,  such  as 


FIG.  9ii  Hliisllcr  nuidc  growth  rates  for  t\\rd  anisotropy,  ij  = 
and  different  values  of  the  overage  hot-eiectron  energy.  Ihe  plasma 
parameters  are  specified  in  the  te.xi. 


FIG  Wi  Hhntler  mode  gionth  rales  fa  fned  loeiage  hot 
electron  energv,  r  ^  =  1.  and  dtffttent  values  of  the  temperattor 
anisotropx 
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unstable  whistler  waves,  to  provide  a  self-consistent 
picture  of  the  hot-electron  distribution  in  energy  and 
the  RF  power  density  required  to  sustain  it.  This 
property  of  the  present  Fokker-Planck  model  is  in 
contrast  to  the  moment  approach  of  Hamasaki  et  al. 
[11],  as  one  might  expect  from  the  non-Maxwellian 
character  of  the  hot-electron  equilibrium  distribution. 
An  interesting  possibility  is  suggested  by  the.se  results, 
namely  that  a  different  class  of  test  distribution  func¬ 
tions,  similar  to  those  obtained  here,  could  be  used  in 
the  moment  approach  to  obtain  useful  analytic  descrip¬ 
tions  of  ECH  equilibria. 

The  numerical  description  obtained  in  the  present 
work  also  contributes  to  clarifying  the  role  of  non- 
adiabatic  scattering  in  limiting  the  hot-electron  energy 
and  it  reveals  in  quantitative  terms  the  importance  of 
the  microwave  integrity  of  the  cavity  in  which  the 
ECH  is  carried  out.  For  single-frequency  heating  it  is 
apparently  necessary  to  have  a  high  enough  cavity 
quality  factor  to  permit  the  ordinary  mode  field 
strength  to  exceed  by  a  large  factor  the  extraordinary 
mode  electric  field  strength.  Similarly.  UORH  can 
only  be  expected  to  yield  the  high  energy  densities 
observed  in  ELMO  if  the  UORH  field  strengths  can 
approach  values  of  the  order  of  10’  V  cm"',  while  the 
corresponding  absorbed  power  densities  remain  around 
1  W  cm or  less. 

The.se  results  are  of  practical  importance  in  designing 
experimental  ECH  approaches  to  achieving  relativistic- 
electron  plasmas;  but  they  also  suggest  that  the 
reduced  Fokker-Planck  model  can  provide  a  useful 
simplified  description  of  other  plasmas  in  which  ECH 
may  occur,  such  as  magnetospheric  plasmas. 

For  our  immediate  purposes,  the  reduced  Fokker- 
Planck  model  provides  a  detailed,  .self-consistent 
description  of  hot-electron  plasmas  that  can  be  analysed 
for  stability.  In  keeping  with  a  large  body  of  empirical 
observations,  we  find  that  whistler  waves  are  likely  to 
be  unstable  unless  the  temperature  anisotropy  is 
reduced  and  the  hot-electron  energy  is  increased. 
Cyclotron  maser  modes  are  shown  to  be  suppressed  by 
the  cold  background  electrons  that  normally  make  up 
the  majority  species  in  an  ECH  plasma. 

In  a  recent  PhD  thesis.  Garner  [29]  has  studied  the 
characteristics  of  microinstabilities  in  an  ECH  plasma 
produced  in  a  magnetic  well  configuration  by  single¬ 
frequency  heating.  Garner  finds  that  the  hot-electron 
component,  whose  temperature  is  4(X)  keV,  is  remark¬ 
ably  stable;  the  observed  unstable  whistler  waves  are 
apparently  driven  by  a  distinct  group  of  warm  electrons 
at  a  temperature  around  2  keV.  The  third  electron 
component  is  a  cold  population  at  a  temperature 


FIG.  10a.  HTiisrler  made  lemparal  grow  th  rates  for  a  Fokker- 
Planck  ECH  equilibrium  w  ith  ad  hoc  values  of  the  temperature 
anisotropy. 


FIG.  lOh.  Hhistler  mode  spatial  growth  rates  and  propagation 
vector  for  the  Fokker-Planck  ECH  equilibrium  with  r)  =  J. 

around  100  eV.  The  cold-,  warm-  and  hot-electron 
densities  are  roughly  equal  and  are  in  the  range  of 
(1-2)  X  10"  em'V  Garner's  thesis  contains  an 
interesting  speculation  regarding  the  observed  stability 
of  the  hot-electron  group  that  is  not  predicted  by  his 
analyses,  namely  that,  if  stochastic  heating  can  occur 
at  energies  above  the  predicted  threshold  for  super- 
adiabatic  electron  behaviour,  the  resulting  distribution 
of  hot  electrons  might  indeed  be  stable  at  high  temper¬ 
atures.  Our  kinetic  model  of  the  relativistic  healing 
process  assumes  that  stochastic  heating  is  maintained 
by  the  overlap  of  relativistic.  Doppler-shifted 
resonances,  and  our  analysis  of  the  resulting  equilibria 
demonstrates  the  relativistic  siabili/ation  of  whistler 
waves  as  well  as  the  cold-electron  stabilization  of  the 
fast-wave  branch  of  the  electromagnetic  mode.  The 
distinct  warm-electron  component  can  be  minimized 
or  even  eliminated  b\  tlie  use  fd  O-niodc  heating  or 
UORH.  as  shown  by  the  present  Fokker-Planck 
results. 
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APPENDIX 

In  brief,  our  picture  of  particle  balance  in  mirror- 
conPined  ECH  plasmas  is  as  follows: 

A-1.  Low-energy  processes 

At  energies  «  S  100  eV,  particle  balanct  is  main¬ 
tained  by  three  dominant  processes,  namely  (i)  ioniza¬ 
tion  of  ambient  neutral  gas,  (ii)  free  flow  of  cold  ions 
along  the  magnetic  lines  of  force  and  out  the  ends  of 
the  magnetic  mirror,  and  (iii)  electrostatic  confinement 
of  electrons  in  an  ambipolar  potential  well  whose 
depth  is  governed  by  the  electron  heating  rate  e  and 
the  cold-ion  escape  rate. 

Equating  the  ion  and  electron  confinement  times 
yields 

e$  =  Lt/Vj 

where  Vj  is  the  cold-ion  flow  velocity  and  L  is  the 
length  of  the  plasma  column.  Furthermore,  the  escape 
of  these  electrons  and  cold  ions  is  balanced  by  ioniza¬ 
tion,  so  that 

Vj/L  *  n()(  U,oni7,ilion^€' 


For  typical  values  of  the  ECH  electric  field. 

Erf  *  10-20  V-cm'',  the  central  magnetic  field. 

B,,  =  3  kG.  the  electron  density ,  n,.  =  5  x  10"  cm 
the  overall  mirror  ratio,  M  =  2.  and  the  mirror  ratio 
at  the  test  electron  turning  point.  M,  =  1.7,  the  criti¬ 
cal  electron  energy  for  runaway,  e,.  can  be  just  above 
100  eV.  (Here.  K  is  the  elliptic  integral  associated 
with  the  electron  bounce  motion.) 

Simultaneously,  with  ECH.  small-angle  Coulomb 
collisions  scatter  heated  electrons  into  the  loss  cone  at 
a  rate  given  by 


(^ )  =  ■ 


5  X  10-* 


n^  In  A  eV- -  cm’ 


The  corresponding  rate  of  pitch-angle  diffusion 
associated  with  ECH  is  given  by 


TT  eElf 


8  BpC 


\  M  I  / 


where  no  is  the  density  of  neutral  gas  in  the  plasma. 
Under  typical  conditions,  the  characteristic  rates  for 
these  three  processes  are  around  5  x  10^  s"'. 

A-2.  Intermediate-energy  processess 

In  the  energy  range  100  eV  s  e  £  10  keV,  the 
dominant  dynamical  processes  are  Coulomb  scattering 
and  ECH.  The  Coulomb  scattering  cau.ses  both 
slowing-down  and  pitch-angle  scattering  of  heated 
electrons.  For  energy  balance,  the  slowing-down  must 
be  offset  by  ECH  for  energies  above  a  relatively  low 
critical  energy  «  2  e|  *  Be^*: 

|rRF(*l)l  =  Ucoul(*l)l 

Using  rudimentary  formulas  for  each  rate,  we  obtain 


where  M,„  is  the  mirror  ratio  at  the  resonant  surface. 
Since  *  f”'-  "hile  Ocnui  “  f RF-driven  pitch- 
angle  diffusion  may  dominate  at  high  energies,  f  >  r j 
=  10  keV,  typically.  In  the  intermediate  energy  range, 
€|  <  c  S  {2.  Coulomb  scattering  is  typically  the 
dominant  loss  process.  The  resulting  pitch-angle 
displacement  of  an  electron  with  initial  pitch  angle  0, 
can  be  estimated  from 


Since  crf  is  almost  independent  of  energy  in  this 
range,  we  have 
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Typically,  <  1,  provided  the  heating  rate  is 

sufficiently  large,  for  example  €rf  <  10^  eV  s"'. 

Some  of  the  runaway  electrons  will  indeed  be  lost, 
especially  for  e  2  «i,  but  the  fractional  hot-electron 
density  obtained  from  energy  balance  alone  provides  a 
useful  estimate  of  an  upper  bound  on  this  ratio.  Note 
also  that  in  the  one-dimensional  Fokker-Planck  model 
we  define  the  hot-electron  density  as  consisting  only  of 
electrons  with  energies  above  a  specified  value,  which 
is  usually  chosen  as  20  keV.  For  e  >  20  keV,  the  RF- 
driven  pitch-angle  diffusion  typically  exceeds  pitch- 
angle  scattering  due  to  Coulomb  collisions;  the  loss 
rate  of  heated  electrons  in  this  energy  range  is  esti¬ 
mated  to  be  very  small.  Since  we  specify  the  total 
electron  density  in  the  one-dimensional  model  and 
restrict  the  hot-electron  component  to  those  electrons 
for  which  Coulomb  scattering  is  negligible,  we  expect 
the  relative  hot-electron  density  obtained  from  energy 
balance  considerations  to  be  a  valid  upper  bound  and  a 
close  estimate. 

A-3,  High-energy  processes 

In  the  energy  range  «  >  10  keV,  Coulomb  scattering 
pitch-angle  diffusion  into  the  loss  cone  can  typically  be 
offset  by  RF-driven  pitch-angle  diffusion  in  which 
heated  electrons  diffuse  towards  the  resonant  surface 
(where  wrf  =  NeB/>m,  N  =  I,  2.  3,...).  Heating  at 
the  second  and  higher  harmonics  of  the  relativistic- 
electron  gyrofrequency  plays  a  major  role  in  the  main¬ 
tenance  of  efficient  stochastic  heating  and  causes 
further  pitch-angle  diffusion  towards  the  midplane  of 
the  magnetic  mirror,  d  =  it  12.  Particle  losses  in  this 
energy  region  are  very  small,  and  energy  balance  is 
dominated  by  relativistic  ECH  and  synchrotron  radia¬ 
tion  up  to  the  energy  at  which  electrons  cease  to  be 
adiabatic.  The  escape  of  a  very  small  number  of  very 
high  energy  electrons  has  little  effect  on  particle 
balance,  but  can  completely  offset  continued  heating  to 
higher  energies. 

In  summary,  we  suggest  that  the  justification  of 
separating  energy  balance  and  particle  balance  condi¬ 
tions  rests  on  the  smallness  of  (A0/0),  as  estimated 
above,  as  well  as  on  the  unusual  nature  of  the  RF- 
driven  pitch-angle  diffusion  that  tends  to  drive  heated 
electrons  back  towards  the  resonant  surface.  We  have 
attempted  to  minimize  the  impact  of  particle  balance 
processes  on  our  one-dimensional  Fokker-Planck 
model  by  regarding  only  electrons  with  energies  above 
20  keV  as  ‘hot  electrons’  and  .specifying  the  total 
density  ab  initio. 
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SECTION  3 


Amplification  of  whistler  waves  propagating  through  inhomogeneous, 
anisotropic,  mirror-confined  hot-electron  plasmas 

G.  E.  Guest  and  R.  L.  Miller 

Applied  Microwavi'  Plasma  Coiiccpis.  hic.,  Carlsbad,  California  92009 
(Received  15  April  1988;  accepted  18  August  1988) 

A  fully  relativistic  local  dispersion  relation  for  whistler  waves  has  been  solved  at  closely  spaced 
points  along  the  magnetic  field  lines  of  a  2: 1  magnetic  mirror  in  which  a  highly  anisotropic, 
spatially  inhomogeneous,  hot-electron  plasma  is  confined.  The  limiting  plasma  parameters  for 
convective  (spatial)  growth  have  been  determined  numerically  and  used  to  identify  plasma 
conditions  leading  to  maximum  amplification  of  input  microwave  signals  introduced  in  the 
form  of  whistler  waves.  The  maximum  gain  has  been  evaluated  numerically  for  a  range  of 
values  of  the  hot-electron  plasma  within  which  all  major  stability  criteria  are  satisfied.  Very 
high  gains  (  ~40  dB)  are  indicated  over  the  entire  range  of  beta  investigated. 


I.  INTRODUCTION 

It  is  well  known  that  whistler  waves  can  be  unstable  in 
mirror-confined  hot-electron  plasmas  because  of  the  pres¬ 
sure  anisotropy  inherent  in  magnetic-mirror  confinement.' 
Indeed,  the  phenomenon  of  unstable  whistler  waves  has  of¬ 
ten  been  invoked  to  interpret  spontaneous  rf  emissions  from 
laboratory*  and  space  plasmas,'  and  has  recently  been  pro¬ 
posed  as  a  way  of  transforming  the  energy  stored  in  mirror- 
confined  hot-electron  plasmas  into  giant  pulses  of  high- 
powered  microwave  radiation.'*'  The  proposed 
transformation  results  from  the  amplification  of  whistler 
waves  propagating  along  magnetic  lines  of  force  and  passing 
through  the  anisotropic,  liot-cleciron  plasma.  Here  we  in¬ 
vestigate  the  maximum  gain  that  a  whistler  can  undergo  in  a 
spatially  inhomogeneous  plasma,  confined  in  a  magnetic- 
mirror  field  of  finite  spatial  extent. 

Because  the  whistler  instability  mechanism  depends  on 
the  kinetic  properties  of  the  hot-electron  equilibrium,  par¬ 
ticularly  the  distribution  of  electrons  in  velocity  space,  theo¬ 
retical  analyses  of  the  growth  typically  use  idealized,  spatial¬ 
ly  uniform  models  of  the  plasma  in  order  to  treat  the 
velocity-space  issues  realistically.'  Under  these  idealized 
conditions,  growth  of  the  whistler  wave  can  be  either  con¬ 
vective  or  absolute,  depending  on  the  plasma  parameters.'’ 
Absolutely  unstable  or  temporally  growing  waves  are  ex¬ 
pected  to  spread  throughout  the  unstable  plasma  and  grow 
in  amplitude  until  quasilinear  or  nonlinear  processes  balance 
the  destabilizing  conditions.  By  contrast,  the  growth  ofcon- 
vectively  unstable  or  spatially  amplifying  waves  may  be  lim¬ 
ited  by  the  finite  transit  time  of  a  wave  packet  across  the 
finite  spatial  extent  of  the  unstable  plasma  as  well  as  by  the 
spatial  nommiformity  of  the  confining  magnetic  field,  rather 
than  quasilinear  or  nonlinear  mecliamsms. 

In  the  present  study  we  use  linearized  theories  to  deter¬ 
mine  the  overall  gain,  Giai),  experienced  by  a  pulse  of  mi¬ 
crowave  power  launched  at  one  end  of  a  mirror-confined, 
hot-electron  plasma  and  |uopagalmg  along  the  magnetic 
lines  of  force  in  the  form  of  w  histler  waves; 

<7(w)  =  1  (la.z  -  )  L)\/\E —  L)\. 

I  lere  o>  is  the  ( radian )  frequency  of  the  wave  launched  at  the 
axial  position  -  L  with  rf  electric  field  strength  E,„. 


After  a  single  transit  through  the  plasma,  the  wave  arrives  at 
the  axial  position  z  =  Z.  with  amplitude  .  Since  this  pro¬ 
cedure  will  yield  the  desired  gain  only  if  the  unstable  growth 
is  convective,  we  have  carried  out  a  detailed  numerical  eval¬ 
uation  of  the  threshold  conditions  for  absolute  growth  and 
restricted  the  range  of  plasma  parameters  to  lie  entirely  with¬ 
in  the  convectively  unstable  range. 

Apart  from  transients  at  the  beginning  and  end  of  each 
pulse,  the  input  wave  amplitude  is  assumed  to  be  constant  in 
time  throughout  the  pulse;  and  the  spatial  and  temporal  be¬ 
havior  of  the  wave  can  be  determined  using  Fourier  and  La¬ 
place  transformed  equations  in  the  conventional  way.® 
When  inverting  the  transformed  variables  it  is  necessary  to 
displace  the  Laplace  integral  path  to  the  real-a>  axis  while 
deforming  the  Fourier  inversion  path  to  prevent  any  zeros  of 
the  dispersion  relation  from  crossing  the  deformed  path  of 
integration  in  the  complex  k  plane,  where  A:  is  the  (parallel) 
wavenumber.  As  is  well  known,  this  procedure  is  only  possi¬ 
ble  if  the  whistler  wave  growth  is  convective;  in  the  case  of 
absolute  growth,  two  poles  from  opposite  halves  of  the  com¬ 
plex  k  plane  merge. 

To  ensure  that  we  consider  only  convectively  unstable 
cases,  we  have  evaluated  the  threshold  conditions  for  abso¬ 
lute  (temporal)  growth  for  the  hot-electron  equilibria  under 
investigation.  For  hot-electron  plasmas  created  by  electron 
cyclotron  heating  (ECH),  the  convective/absolute  bound¬ 
ary  can  be  conveniently  represented  in  two  dimensions  by 
fixing  some  of  the  plasma  parameters  in  ways  that  reflect 
fundamental  properties  of  the  ECH  process.  Thus,  for  exam¬ 
ple,  the  total  plasma  density  is  usually  limited  by  microwave 
cutoff  so  that  ai^f„  /a)l  1.  Here  and  are  the  electron 
plasma  and  microwave  frequencies,  respectively.  Moreover, 
the  relative  density  of  hot  electrons  n„/n  =  A  may  be  limited 
by  the  onset  of  flutelike  instabilities  in  simple  magnetic-mir¬ 
ror  configurations;  although  no  such  limitation  occurs  in 
magnetic  "well”  configurations  with  favorable  curvature  of 
the  magnetic  lines  of  force.  Even  in  such  a  flute-stable  config¬ 
uration,  however,  the  plasma  pressure  may  be  limited  by 
equilibrium  conditions  such  as  the  so-called  "mirror  mode,” 
and  the  value  of  hot-electron  beta,  may  be 

limited  in  inverse  proportion  to  the  anisotropy.’  Here  B  is 
the  magnetic  intensity  and  /i„  is  the  permeability  of  free 
space. 
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If  no  other  stability  criteria  are  violated,  the  gain  G(a>) 
can  be  evaluated  up  to  the  threshold  conditions  defined  by 
the  C /A  boundary  in  a  straightforward  way;  this  is  the  main 
objective  of  the  work  described  here.  The  gain  is  then  given 
directly  by 

G(<y)  =  exp^  —  2  J  ki({o,z)dz 

where  k,  is  the  imaginary  part  of  the  wave  vector  corre¬ 
sponding  to  solutions  of  the  local  dispersion  relation  with 
real  frequency  lo. 

In  order  to  account  for  the  z-dependent  equilibria  con¬ 
fined  in  finite  sized  magnetic-mirror  configurations  we 
adopt  simple  model  equilibria  with  hot-electron  distribution 
functions  mapped  along  magnetic  field  lines  through  adiaba¬ 
tic  invariants.  The  properties  of  these  equilibria  are  summar¬ 
ized  briefly  in  Sec.  II.  A  fully  relativistic  dispersion  relation 
is  then  used  to  evaluate  the  C /A  boundary  for  whistler 
waves.  The  dispersion  relation,  while  fully  relativistic,  can 
give  only  a  local  evaluation  of  the  wave  properties,  since  it  is 
derived  under  the  assumptions  of  spatial  uniformity.  This 
procedure  is  expected  to  yield  a  valid  description  provided 
the  scale  length  for  variation  of  the  plasma  equilibrium  is 
much  greater  than  the  wavelength  of  the  whistler  waves. 
Since  we  are  primarily  interested  in  quasioptical  systems 
with  scale  lengths  much  larger  than  wavelengths,  this  condi¬ 
tion  is  not  particularly  restrictive.  In  effect,  we  are  construct¬ 
ing  a  WKB  solution  from  numerical  solutions  of  the  local 
dispersion  relation  within  plasma  parameter  ranges  for 
which  the  unstable  growth  is  guaranteed  to  be  convective. 
The  numerical  analysis  of  the  dispersion  relation  and  typical 
results  are  recapitulated  in  Sec.  Ill  The  gain  is  then  evaluat¬ 
ed  from  real-rj  solutions  of  the  dispersion  relation  in  a  way 
that  is  discussed  at  some  length  in  Sec.  IV.  Several  cases  are 
presented  there  in  which  the  temperature  anisotropy  and 
hol-electron  beta  are  varied  so  as  to  move  along  a  particular 
C /A  boundary  (corresponding  to  a  fixed  relative  hot-elec¬ 
tron  concentration,  A  =  30%).  Substantial  values  of  the 
gain  are  indicated  for  hot-electron  plasmas  of  typical  labora¬ 
tory  dimensions. 

The  major  findings  arc  summarized  in  Sec.  V,  and  the 
nature  of  these  findings  is  discussed  In  particular,  the  unre¬ 
solved  issue  of  the  role  of  quasilinear  or  nonlinear  processes 
is  given  some  emphasis  because  of  the  indications  that  the 


maximum  gain  predicted  in  the  line.ar  llienry  varies  onlv 
weakly  with  beta  and  remains  l.irpe  in  plasmas  with  high 
energy  densities,  where  satmation  levels  of  output  power 
might  be  expected  to  be  very  high. 


II.  THE  MODEL  EQUILIBRIA 

In  this  section  we  summarize  the  way  in  which  the  equi¬ 
librium  properties  (vf  the  plasma  depend  on  distance  along 
the  magnetic  lines  of  force.  Local  prc'perlies  of  election  cy¬ 
clotron  heated  hot-electron  plasmas  have  been  discussed  ex¬ 
tensively,^  and  a  few  studies  have  attempted  to  model  the 
axial  dependence  in  nonrelativistic  formulations  of  the  hot- 
electron  dynamics.  Here  we  follow  the  conventional  ap¬ 
proach  to  mapping  relativistic-electron  equilibria  along 
magnetic  lines  of  force  assuming  the  inv  ariance  of  the  total 
energy  f,  and  the  magnetic  moment  //: 

f  —  mc'(  y  —  I )  --  rnc  [  I  i  ( u  /c  )  -  I  J 
and 

fi  =  [mil]  )/{2B). 

In  these  expression.  is  the  relativistic  I  orentz  factor  and  u 
is  the  momentum  per  unit  rest  mass  with  components  u  and 
u  perpendicular  and  parallel,  respectively,  lo  the  magnetic 
field  B.  Here  m  is  the  electron  rest  mass  and  c  is  the  speed  of 
light  in  vacuum. 

The  relative  intensity  of  the  sle.'idy  magnetic  field  is  as¬ 
sumed  lo  vary  along  the  line  of  force  of  interest  as  follows; 


Here  z  =  0  defines  the  so-called  “midplane"  on  w  hich  the 
magnetic  intensity  reaches  its  minimum  value  />  —  1.  while 
z  —  Z.  defines  the  so-called  "mirror  throat,"  where  5  reaches 
its  maximum  value  h  =  M.  the  mirror  ratio 

The  cold-electron  component,  which  typically  cont.iins 
the  majority  of  elections,  will  he  modeled  as  an  isotropic, 
pressurelcss  fluid  w  hose  density  is  independent  of  z  By  con¬ 
trast.  the  hol-electron  component  is  assumed  lo  be  aniso¬ 
tropic  and  relativistic  and  described  by  a  distribution  func¬ 
tion  of  the  form" 


If  u„  is  the  momentum  per  unit  rest  mass  evaluated  at  the 
midplane  z  =  0,  then  in  terms  of  the  constants  of  motion, 

wi')ii  =  (z)  -f  (A  -  1 )  [  uj  (z)  ]/h 

(2) 

and 


u;„  -  2/in„/m  -  u\/h.  (3) 

Merc  1/,,^  and  i;,,.  are  the  compi'ncnls  of  u,,  parallel  and  per¬ 
pendicular.  respectively  . !('  the  static  magnetic  fielil 

A  z-dependeni  loss  cone  in  u  space  is  defined  by  setting 
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n^(j|/r>^(OI 


(4) 


=0: 


ujj  ( Z. )  =  u’  -  /m  =  0 

=  Mj5  (2)  +  u]  (z)  -  t/f  (z)B,.,,yB(z), 


or 


^  _  A/ 
5(^)  6 


-  I. 


If  we  now  collect  terms  in  t/J  (z)  in  the  exponential  factor  we 
have 

=  2  Ci'-(Ayexp(  -  4  - 

/-n  \ba^  J  \  n,,  a\(z)J 


if  — 1-< - , 

h  -  1 


(5) 


=  0.  if^>-  ^ 


where 

of  (2)  =  (a;„a||h)/{a;„(h  -  1)  +  (tj  |, 
or,  if  ?7  =  a;n/ai;, 

af  (2)  =  r)bal/[-q(h  -  I )  +  I  ]. 

The  2-dependent  hot-electron  density  is  now  given  by 


n^z) 


(6) 


Jr-u'  ( *(  M  \  ha\  )  \  a\  (z)J 

and  for  the  simple  bi-Gaussian  distribution,  /  =  0,  we  have 


n^(z)  =  ir'''a\  (z)a,  €!"'{  I  +  (of/or'  (2) ) 
X\b/(M -b)]} 
or 


(8) 


FIG  I  Varianon  of  relative  hot  electron  density  and  magnetic  inlensily.  h. 
with  distance  along  a  magnetic  field  line 


w^(2)  (2)/  ^  {ni/a]„)  -  D)  y- 

Om,  \  1  -i  [o'-Zo-;  (z)  ]  |/)/(  U  h)\) 

, - 1 -  [J^  m 

I)  t  i-y  (M  !)• 

An  illustrative  case  is  show  n  in  I  ig  I  lot  a  tnit lor  ratio 
Af  —  2  and  two  ditTerent  tempeialurc  anisotropies,  t/  —  5 
and  20.  The  2-depetidcnt  relative  magnetic  inicttsity  b  is 
shown,  together  with  the  z-dcpendetit  relative  hot-clectroti 
densities  for  the  two  ditTerent  tetnpertiture  anisolropi  .'s.  I  he 
figure  indicates  the  degree  to  w  hich  higher  temperature  atii- 
sotropy  is  associated  with  increased  peakittg  of  the  hot-elec¬ 
tron  density  at  the  midplane.  This  property  of  the  equilibria 
will  be  shown  to  have  a  major  impact  on  the  overall  gain  of  a 
whistler  wave  propagating  along  the  magnetic  line  of  force. 
In  quantitative  terms,  the  z  dependence  of  hot-electron  den¬ 
sity  is  given  by 

d  _  db  (  ( I  -  ?;)  I  {M  -  b) 

dz  nJO)  ~  d2  V|l  -I-  7;(/>  -  /)p  yj  T.V/  -  1 ) 

— — i - J - ) 

2v(Af-  l)(Af-  h)  1  Vib  -  I)/ 

db/ ,  I  \ 

-  — I  I  -  V  -  --  - -  .  as  h-.  1  ( 10) 

dz\  2(Af  I)/ 

Although  it  is  tiot  incluiled  in  l  ig  I,  the  Ictnpcraturc  att- 
isotropy  also  decreases  with  distattcc  from  the  midplanc  as 

i)b 

n,;  1  4  >/(/’  ■  I ) 

so  that 

da]  _  i}b'(]  -  I/) 

dz  1 1  +  v(^  -  I )  I'  *  "  * 

Thus  both  temperature  anisotropy  and  the  relative  hot-elec- 
tron  density  decrease  with  distance  away  from  the  midplane, 
so  that  maximum  amplification  will  occur  at  the  midplane. 

In  contrast  to  this  strivng  z  dependence  of  the  hi't-elec- 
tron  density,  the  cold-electron  densitv  is  expected  to  remain 
nearly  constant  in  z,  except  for  the  spatial  compression  asso¬ 
ciated  with  the  increasing  magnetic  intensity.  The  real  value 
of  k((o)  is  determined  mainly  by  the  cold-plasma  density 
and  the  magnetic  intensity,  and  is  slowly  varying  along  the 
magnetic  field  lines. 


III.  CONDITIONS  FOR  MAXIMUM  GAIN 

The  maximum  amplification  of  an  externally  launched 
whistler  wave  is  likely  lobe  acliie-.  cd  in  hot-electron  plasmas 
that  are  just  below  the  threshold  for  absolute  or  temporal 
growth.  This  threshold  is  frequently  referred  to  as  the  C /.I 
boundary,  which  is  a  hypersurface  in  a  space  spanned  by  the 
plasma  parameters  On  one  side  of  this  houndarv  onlv  con¬ 
vective  or  spatial  growth  is  possible,  vvhilc  c'li  the  I'pposiic 
side  of  the  boundary,  absolute  or  temporal  grow  th  occurs 
Thus,  in  order  to  exphne  the  maximum  gain  that  can  be 
achieved,  it  is  first  necessary  to  define  the  ( '  /,  I  boundary  in 
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terms  of  equilibrium  parameters  representative  of  physically 
realizable  plasmas.  This  section  summarizes  the  results  of  a 
numerical  determination  of  the  C  /A  boundary  for  unstable 
whistler  waves  in  hot-electron  plasmas  whose  local  proper¬ 
ties  are  modeled  by  the  class  of  infinite,  homogeneous,  rela¬ 
tivistic-electron  distribution  functions  cited  earlier,  namely. 


The  general  methodology  adopted  here  for  determining 
the  C /A  boundary  has  been  discussed  at  length  by  Briggs'" 
and  consists  ultimately  in  locating  simultaneous  solutions  of 
the  dispersion  relation  and  its  derivative  with  respect  to  k, 
the  (complex )  wavenumber,  for  real  values  of  the  frequency 
<u.  This  general  methodology  is  implemented  in  our  case 
through  novel  numerical  techniques  appropriate  to  the  class 
of  distribution  functions  used  here.  These  techniques  have 
been  fully  described  by  Guest  and  Miller.'’ 

It  has  been  well  established  by  many  investigators'  that 
the  grow’th  rate  of  unstable  whistler  waves  increases  with 
increasing  values  of  the  temperature  anisotropy  and  the  hot- 
electron  beta,  where  beta  is  the  ratio  of  the  kinetic  pressure 
nT  to  the  magnetostatic  pressure  B'/l/i,,.  For  hoi-eleciron 
plasmas  created  by  electron  cyclotron  heating  ( ECH ),  sev¬ 
eral  fundamental  aspects  of  the  heating  and  confinement 
processes  limit  the  physically  accessible  values  of  beta  and 
anisotropy  in  ways  that  can  usefully  restrict  the  range  of 
plasma  parameters  to  be  considered  here.’  ’  For  example, 
effective  penetration  of  the  ECH  power  typically  limits  the 
maximum  (total)  electron  density  to  values  such  that 
0)1^  where  is  theangular  frequency  of  the  ECH  pow¬ 
er.  Typically,  the  relati  e  fraction  of  electrons  heated  to  high 
energies,  A  —  «*/«,  is  much  less  than  unity.  Idealized  theo¬ 
ries  of  flutelike  instabilities  have  been  used  to  estimate  criti¬ 
cal  values  of  A  set  by  thresholds  for  occurrence  of  flutelike  or 
ballooning  modes:  A^,„  —0.3  Relativistic  cflects  appear  to 
increase  this  substantially,  but  we  will  usually  assume 
A<0.5.  In  small  magnetic  mirrors,  nonadiabatic  effects  can 
limit  the  maximum  hot-electron  energy  to  values  such  that 
/7<0.05F.  Here^is  the  (relativistic)  electron  gyroradiusand 
Z.  is  a  characteristic  scale  length  of  the  magnetic  geometry. 
The  condition  for  equilibrium,  sometimes  referred  to  as  the 
mirror  instability  criterion,  limits  the  product  of  anisotropy 
and  beta,  roughly  as  )9( »/  —  /)  <  I .  For  the  anisotropic,  hot- 
electron  distribution  functions  used  in  the  present  work,  a 
more  accurate  estimate  of  the  beta  limit  corresponding  to  the 
mirror  criterion  can  be  obtained  from  the  more  general  con¬ 
dition  for  existence  of  equilibria'": 


B  dB 


In  Ihenonrelativi'^ficapproximation./i, ,  the  component 
of  plasma  pressure  perpendicular  to  the  static  magnetic  field, 
is  given  by 


=  mn(vl)  =  mn(vl)  =  mti(r]/2) 

=  TT J  J  I',  <fr,/(i',.r||)wr|. 


Herei'’,,,  =  r'h/iM  -  />)  deliiics  the  loss  cone  at  the  point  r 
where  the  local  magnetic  intensity  is  /)  -  />/?,,  Note  that 
I’l/n  — Oas  .3/—  cc.  If  we  use  the  bi-Maxwellian  distribution 
function  corresponding  tc'  /  =-  O"  and  evaluate  the  derivative 
ofp,  at  the  midplane,  we  obtain  from  the  mirror  criterion 

^  ^ _ ' _ 

(M  -  i)(v  -  i){  ’  (  M  I  )  hi  '/H  ' 

(  14) 

where  r]  =  a^,/a^^.  We  have  evaluated  the  limiting  value  of 
/?,  ( r})  for  A/  —  2  and  displayed  this  boundary  in  I'ig.  2.  Note 
that  as  A/—  oo  this  criterion  approaches 

ffiiV  -  I) <  1 . 

In  Fig.  2  we  show  the  value  of  anisotropy  on  the  convec¬ 
tive/absolute  boundary  r;,  ^  as  a  function  of/?,  for  condi¬ 
tions  that  are  appropriate  to  an  ECH-generated  hot-electron 
plasma  We  have  set  w’,  Au’,  =  1  and  A  =  0,3.  and  varied  ;■ 
to  cover  the  indicated  range  of  values  of/9; .  The  associated 
values  of  y  are  shown  by  the  curve  labeled  y  -  1 .  Here  y  is 
the  average  of  y  over  the  equilibrium  distribution  function. 
Eq.  (5 ),  in  which  only  the  /  =  0  component  is  retained.  1  he 
approximate  mirror  instability  criteria  fSiv  —  I )  =  I  and 
r  =  0  are  also  shown  here;  note  that  the  C /A  boundary  lies 
slightly  above  the  equilibrium  limit  for^,  >  \0^r-20^r.  The 
nonrelativistic  C /A  boundary  given  by  liy  oshi  cl  al. ' '  is  also 
shown  for  comparison  with  the  fully  relativistic  numerical 
results.  For  a  given  value  of/?,  the  relativistic  expressions 
indicate  that  slightly  higher  values  of  anisotropy  are  re¬ 
quired  for  absolute  growth 

The  C /A  boundary  shown  itr  Fig.  2  satisfies  the  condi¬ 
tions  for  physical  accessibility  identified  here  for  /?  •  lO'T  - 
20%  and  thus  will  be  used  to  estimate  the  maximum  gain 
that  can  be  achieved  at  lower  values  of  beta.  The  results  of 


f1(i  2  V.ni;ili«»n  «’f  t»f  lln  f  f  htMinvI.UN  >/,  , .  ;»inl  .i\ cf  .tpc 

h<il  clcclr«»n  cncigN.  ) ,  wiih  /?  .  tht  (vm pt-iidK  iiI.ii  fn'i  cln  Mi»ti  N-i.i  llcio 
ri’  /»•»'  -  1  jnd  A  (1  '  AIm>  sfunMi  I*,  llu*  niiMtH  tliicsht'ld 
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thi’  *v  ‘!uation  of  gain  are  described  in  Sec  IV, 

vVe conclude  the  present  section  with  one  furlhcr  set  of 
results  on  the  local  wave  properties  of  unstable  whistlers  at 
the  C /A  boundary  shown  in  Fig.  2.  These  are  displayed  as 
functions  of  in  Fig.  3,  where  we  see  that  both  the  real  and 
imaginary  components  of  the  propagation  vector,  together 
with  the  wave  frequency,  are  decreasing  functions  of  beta. 
This  decrease  reflects  two  underlying  facets  of  the  C/A 
boundary  along  which  the  parameters  are  evaluated,  name¬ 
ly,  that  the  maximum  anisotropy  for  convective  growth  is 
decreasing  with  increasing  values  of  beta,  and  that,  in  the 
present  case,  beta  is  being  increased  by  increasing  the  aver¬ 
age  electron  energy  and  hence  the  relativistic  mass  shift.  The 
relativistic  effects  can  be  removed  by  holding  the  hot-elec¬ 
tron  energy  fixed  and  varying  beta  by  varying  A.  This  proce¬ 
dure  has  been  followed  in  generating  the  results  shown  in 
Figs.  4  and  5. 

Figure  4  again  shows  the  anisotropy  at  the  C /A  bound¬ 
ary  as  a  function  of  beta,  together  with  the  corresponding 
values  of  A  and  the  mirror  instability  criterion  for  y  =  1.4. 
This  value  of  y  corresponds  to  an  average  hot -electron  ener¬ 
gy  just  above  1(X)  keV  (T  ~  111  keV).  Note  that  absolute 
growth  is  only  possible  for  /?,  >0  10  corresponding  to 
A  =  0. 1 55,  and  that  the  C /A  boundary  drops  below  the  mir¬ 
ror  stability  limit  for^,  >  0, 17  or  A  >0.27.  The  wave  param¬ 
eters  along  the  C /A  boundary  are  show  n  in  Fig.  5,  and  here 
reflect  only  the  decreasing  value  of  anisotropy,  as  well  as  the 
increasing  value  of  A  In  the  range  that  satisfies  the  mirror 
stability  criterion,  the  wave  parameters  decrease  only  slight¬ 
ly.  suggesting  that  relativistic  effects  are  mainly  responsible 
for  the  large  decrease  in  wave  parameters  shown  in  Fig.  3. 

IV.  EVALUATION  OF  MAXIMUM  GAIN 

We  now  estimate  the  maximum  achievable  gain  for  var¬ 
ious  plasma  parameters  by  numerically  evaluating  the  inte- 


r!G  4  The  C/A  boundary  for  fixed  y  ~  14  and  Ihc  relative  hoi-elcclron 
density.  A.  together  with  the  mirror  instahilit)  criterion 

gral  along  magnetic  field  lines  of  the  imaginary  part  of  the 
propagation  vector: 

lnG(rt;)=  k,{(i).z)dz.  (15) 

This  procedure  is  valid  provided  the  growth  is  convective. 
We  will  examine  a  sequence  of  cases  in  which  the  plasma 
beta  is  increasing,  and  for  each  value  of  beta  w  e  w  ill  assume 
the  maximum  anisotropy  for  which  the  whistler  growth  is 
convective.  Since  the  plasma  density  and  anisotropy  are 
greatest  on  the  midplane,  it  is  the  midplane  values  of  plasma 


Flo  3  Wavt  parameters  versus^  for?;,  ,  (m  the  C/.-thoiindar)  shown  in 
Fig  2 


I'lO  s  Wase  pai.imclcrs  along  Ihc  (  /.r  houiiil.iis  foi  Ihc  cases '.how n  m 

Fig  4 
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and  wave  parameters  on  the  C/A  boundary  that  determine 
the  maximum  anisotropy  and  the  corresponding  wave  fre¬ 
quency  for  each  value  of  beta  The  anisotropy  and  density 
are  readily  determined  at  each  point  along  the  magnetic  field 
line;  and  the  dispersion  relation  is  then  solved  locally  to  ob¬ 
tain  the  complex  propagation  vector  at  that  point  for  the 
frequency  specified  at  the  midplane.  The  gain  integral  is  then 
evaluated  numerically. 

We  have  carried  out  this  evaluation  of  the  gain  for^, 
=z5%,  10%,  15%,  20%,  and  25%  in  the  case  shown  in  Fig. 
2:tuJ,/<uJ,  =  1,  A  =  0.30,  and /  =  0.  The  magnetic  geometry 
assumed  was  a  fixed  M  magnetic  mirror;  and  thus  no 
attempt  was  made  to  relate  the  equilibrium  value  of  the  an¬ 
isotropy  to  the  shape  of  the  magnetic  field.  We  return  to  this 
point  in  Sec.  V.  The  midplane  parameters  for  the  cases  exam¬ 
ined  are  tabulated  in  Table  1.  The  series  was  terminated  at 
=  25%  since  the  mirror  instability  criterion  was  not  satis¬ 
fied  at  the  C /A  boundary  for  greater  values  of  0,. 

The  last  column  of  Table  1  gives  the  value  of 


The  gain  is  then  given  by 

=  exp[4(L(d,,t,/C)I  ]. 

For  typical  conditions,  the  applied  microwave  power  reso¬ 
nates  with  the  second  harmonic  of  the  local  electron  gyrofre- 
quency  at  the  midplane.  Thus 

<Uc.o  = 

and  we  can  estimate  the  scale  factor  as  follows: 

4irLf„  _  4nL 
c  ~  c  ' 

where  is  the  free-space  wavelength  of  the  microwave 
heating  power.  In  order  of  magnitude,  L//1„>C7(!0)  for 
typical  existing  experimental  devices.  The  indicated  values 
of  gain  are  then,  again,  in  order  of  magnitude, 

C{a))~C7(10^)~40dB. 

In  addition  to  the  peak  growth  rates  show  n  in  Table  I  we 
can  readily  estimate  the  effective  bandwidth  by  calculating 
the  total  gain  for  whistler  waves  of  varying  frequencies  prop¬ 
agating  through  a  particular  anisotropic,  mirror-confined, 
hot-electron  equilibrium.  For  definitiveness,  we  again 
choose  the  plasma  parameters  to  be  near  the  convective/ 


MO  6  1  Ilf  frfiiiifiiis-dcpfnflfnt  p:iin  f.iflnr.  l'(i  i"'  ,  fi'r 

whistler  waves  of  varving  frcquencN  m  ht'i  ifeiiron  f(|mfihria  wHh  feni- 
petal urc  anisotropies  near  ihc  coiuci  Ms  p  'ahsi'luic  htMiiul.irv  aiul  /?  . 

.  and  20‘T-  Plasma  and  wave  paianuMcrs  ft'r  these  three  t  .ises  are  IrstecI 
in  Table  I. 


absolute  boundary  for  the  cenlci  ficqiiency  ('f  the  unstable 
frequency  band.  lllustrati\c  cases  of/?,  f  5%,  10%,  and 
20%  are  shown  in  Fig.  6,  where  the  dimensionless  gain  fac¬ 
tor. 


is  plotted  against  thefiequency  noiinalixed  to  the  center  fre¬ 
quency  Note  that  the  bandwidth  increases  as/?  increases, 
i.e.,  as  the  average  hot-electron  energy  increases  I  he  full- 
width  at  half-maximum  varies  from  roughly  0  15  at  0 
=  5%  to  more  than  0.25  at  /?  =  20%. 

Thus  the  linear  theory  suggests  that  very  high  gain  is 
possible.  It  should  be  emphasized  that  the  linear  theory  can- 


TABLE  I.  Summary  of  Ihc  wave  and  plasma  parameters  .u  ihc  C /A  Nnindary  for  varying  values  of  /?  l  or  all  these  vases,  r.r.  /r  >  I  ami  A  H  V) 


/?.(%) 

V,  /. 

rr/Ai'..,, 

(■'  *  <  * 

.1  1  „  /, 

5 

15  26 

0  699R 

1  623 

0  7267 

1  (l'>1 

0  07(3 

10 

It  83 

0  5769 

1  116 

0  5964 

1  |97 

0  O'  1 

IS 

9  61 

0  4S40 

I  (KM 

0  5t50 

1  1I2V 

on'^0 

20 

7  972 

0  tO'IS 

()qiS7 

0  4487 

1  4'7S 

tl  06 

25 

6  698 

0  3484 

ofiioq 

0  1942 

1  ''2 

0  0(4 

30 

5.724 

0  2977 

0  7097 

0  1485 

1  71' 

35 

4.988 

0  2562 

0  6280 

0  11 0(1 

t  810 

40 

4  426 

ISJ 

0  5617 

0  27K7 

? 
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not  provide  any  indication  of  the  maximum  output  power  or, 
equivalently,  the  saturation  value  of  the  input  power.  These 
must  be  obtained  from  quasi  linear  or  no;)linear  models  of  the 
energy  transformation  process. 

It  is,  in  fact,  most  striking  that  the  gain  decreases  only 
very  slowly  with  increasing  beta;  one  would  expect  that  the 
maximum  output  power  would  be  proportional  to  beta.  As 
we  remarked  earlier,  the  fact  that  the  spatial  growth  rate 
decreases  with  beta  is  in  part  attributable  to  the  .somewhat 
arbitrary  assumption  of  fixed  A.  Beta  is  then  proportional  to 
the  average  energy  of  the  hot  electrons,  and  relativistic  ef¬ 
fects  reduce  the  growth  rates.  Perhaps  more  important  to  the 
weak  beta  dependence  of  the  gain  is  the  fact  that  at  low  beta 
the  C /A  boundary  permits  very  high  anisotropy.  High  an¬ 
isotropy,  in  turn,  leads  to  very  large  spatial  growth  rates  near 
the  midplane.  Although  lower  anisotropy  does  lead  to  less 
peaking  of  the  density  near  the  midplane,  the  growth  is 
greatly  reduced  away  from  the  midplane  by  the  increased 
local  magnetic  field  strength;  is  too  small  to  permit 

large  local  growth  rates. 

V.  CONCLUSIONS 

In  Sec.  IV  we  demonstrated  that  the  amplification  of  an 
externally  launched  whistler  wave  can  exceed  40  dll  over  a 
relatively  wide  range  of  values  of  beta,  5%  <25%,  for 

typical  ECH  plasma  parameters:  =  I,  A  —  n^/n 

=  0.3,  and  1.1*.  y<  1 .6.  The  corresponding  limititig  values  of 
temperature  anisotropy  were  15^?;,  ,,  >6.  These  values  of 
ani.sotropy  are  much  larger  than  expected  in  a  typical  ECU 
plasma  confined  in  a  2:1  magnetic  mirror.  Daiull^  has  pro¬ 
posed  the  use  of  adiabatic  compressiiui  to  increase  »;  to  the 
C /A  boundary  value.  We  have  not  attempted  to  include 
adiabatic  compression  in  the  present  studies.  1  he  gain  calcu¬ 
lation  included  the  strong  spatial  variation  in  hot-electron 
density  and  anisotropy,  together  with  the  spatial  variation  in 
the  underlying  2  I  magnetic-mirror  field  We  conclude  from 
these  findings  that  strong  amplification  is  possible  for  repre¬ 
sentative  hot -electron  equilibria  in  the  linear  regime;  i.e  ,  for 
input  signals  below  some  power  level  at  which  quasilinear  or 
nonlinear  effects  will  saturate  the  growth  of  the  whistler 
wave.  Clearly  the  linear  theory  cannot  provide  an  estimate  of 
this  critical  parameter;  some  type  of  quasilinear  calculation 


of  the  saturation  of  whistler  growth  is  required  t<i  address 
this  issue. 

Existing  quasilinc!ir  studies  of  unstable  whistler  w;i\es 
in  infinite  homogeneous  plasmas  have  indicated  that,  within 
that  highly  idealized  model,  roughly  W%  of  the  hot-electron 
energy  can  be  transformed  into  wave  enetgv  tiefoie  the 
growth  rate  is  reduced  to  rero.  If  these  estimates  provide  a 
valid  indication  for  finite,  inhoniogeueous,  mirror-confined 
plasmas,  it  would  appear  entirely  feasible,  by  operating  at 
higher  values  of  beta  in  high  field  strength  magnetic  fields  to 
achieve  very  high  output  power  levels  In  fact,  it  may  well  be 
that  the  infinite  homogeneous  quasilinear  model  prov  ides  an 
underestimate  of  the  saturation  power,  since  resonant  elec¬ 
trons  can  escape  a  real  magnetic  mirror  after  giving  a  sub¬ 
stantial  fraction  of  their  perpendicular  kinetic  energy  to  the 
wave.  In  effect,  a  nearly  constant  level  of  anisotropy  is  auto¬ 
matically  maintained  in  an  actual  magnetic  mirror,  whereas 
the  infinite,  homogeneous  quasilinear  model  allows  the  ani¬ 
sotropy  to  decrease  as  a  result  of  the  wave  growth. 
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4.  RECENT  RESULTS 


In  this  secliori  oje  bunimarize  the  mobl  recent  results  achieved, 
including  a  preliminary  estimate  of  the  maximum  fraction  of  the  plasma 
energy  that  can  be  transformed  into  RF  energy,  and  a  brief  discussion  of 
the  propagation  characteristics  of  uuhistlers  in  candidate  PEMS  devices. 
UUork  in  both  of  these  areas  ujas  undertaken  in  anticipation  of  intensive 
research  into  the  non-linear  limits  on  the  PEMS  concept  and  the 
geometrical  optics  issues  that  must  be  resolved  to  ensure  efficient 
coupling  of  the  output  microiuave  fields. 

THERMODYNAMIC  ESTIMATES  OF  THE  MAXIMUM  FRACTION  OF  THE 
HOT-ELECTRON  ENERGY  THAT  CAN  BE  TRANSFORMED  INTO 

MICROUJAVES 


Thermodynamic  techniques  developed  by  C.  S  Gardner  (6],  T.  K. 
Fouuler  (5)  and  others  can  be  used  to  estimate  the  maximum  fraction  of 
the  mirror-confined  hot-electron  kinetic  energy  that  can  be  transformed 
into  microujeves.  Here  uue  recapitulate  the  essential  elements  of  the 
thermodynamic  approach  and  apply  it  to  the  PEMS  concept. 

The  microuuave  energy  within  some  volume  V  is  given  by 

F  =  Jd^x  (£oE2/2  +  B2/2po) 

V 

=  Jd^x  EqE^ 

V 

Following  Fowler  [5]  we  define  a  free  energy  for  the  anisotropic. 
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hot-eLectron  pLasma: 


ujhere 


ft  H  l/Zjd^x  JcPv  {G(f)  -  6{g}  +  mv^ff  -  g)/2)  +  F 
V 


G(h)  £  hLn(h/C]-h 


and  the  reference  distribution  function,  g,  is  given  by 


9  £  C  exp 


2T 


corresponding  to  thermal  equilibrium.  The  parameters  C  and  T  are 
chosen  to  minimize  the  value  of  ft(0)  at  an  initial  time  uuhen  the 
microojave  energy  is  negligible.  It  can  be  shouun  ( )  that  ft(0)  is  an  upper 
bound  on  F;  and  in  particular,  if  the  equilibrium  hot-electron  distribution 
function  f  is  the  1  =  0  bi-Maxuuellian  case  assumed  in  most  of  the 
studies  reported  here, 


F 


2 


Thus,  averaged  over  the  plasma  volume  V, 


F/V  •  < 


2 


The  total  kinetic  energy  stored  in  the  mirror-confined  hot  electron  is 
given  for  this  distribution  function  by 


UJ  =  nTi 


2yJ 
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Thus,  the  maximum  fraction  of  the  total  hot -electron  kinetic  energy  that 
can  be  transformed  into  microu/ave  energy  is  given  by 

UL  ^  -L  lM-lP 

U)  4y  (y  ♦  0.5). 

In  the  limit  y  oo,  this  fraction  approaches  25%;  the  limit  is  20%  for  y  = 
11.6  and  10%  for  y  =  3.1.  Thus,  over  the  range  of  values  of  temperature 
anisotropy  examined  in  this  report,  10  -  20%  of  the  hot-electron  energy  is 
expected  to  be  transformed  into  microiuave  energy  in  an  individual 
microujave  pulse.  It  is  important  to  recognize  that  the  energy  remaining 
in  the  hot-electron  plasma  is  not  lost,  and  can  be  utilized  in  succeeding 
pulses.  In  fact,  the  repetition  rate  of  the  pulses  can  be  much  greater 
than  the  inverse  of  the  initial  buildup  time  of  the  hot-electron  plasma 
precisely  because  of  the  energy  remaining  after  each  pulse;  it  is  only  the 
energy  transformed  into  microoiave  pouuer  that  must  be  replaced.  The 
average  conversion  efficiency  can  therefore  be  very  high,  provided  the 
rote  of  energy  replocement  (by  UORHJ  exceeds  the  rote  ot  luhich  hot 
electrons  are  lost  from  the  magnetic  mirror  by  Coulomb  scattering  into 
the  loss  cone.  These  losses  can  be  very  small;  the  dominant  pouuer  drain 
is  likely  to  be  that  uuhich  is  associated  uuith  the  continuous  throughput  of 
louu-temperature  plasma.  This  pouuer  is  supplied  through  fundamental 
ECH. 
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GEOMETRICAL  OPTICS  OF  UIHISTLER  (HAVES 
IN  CANDIDATE  PEMS  GEOMETRIES 


UJe  have  carried  out  a  preliminary  study  of  the  propagation  properties 
of  uuhistler  uuaves  Launched  at  one  end  of  a  PEMS  device  and  extracted 
from  the  opposite  end.  Since  the  luave  propagation  is  governed  mainly 
by  the  colder  plasma  component,  uue  can  obtain  a  preliminary  vieuu  of  the 
uuave  characteristics  by  making  ta;o  plausible  assumptions  regarding  the 
variation  of  cold-plasma  density  along  the  magnetic  lines  of  force. 

These  are 

1.  ujithin  the  magnetic  mirror  region  the  cold-plasma  pressure  and 
temperature  are  constant;  and 

2.  in  the  region  outside  the  magnetic  mirror,  the  plasma  flouus  along 
the  diverging  magnetic  field  lines,  and  Local  Ionization  is  negligible. 

The  first  assumption  implies  that  the  dimensionless  ratio  of  electron 
plasma  and  gyrofrequencles,  cope/cojjg,  varies  inversely  luith  magnetic 
intensity,  B; 

COpe/cOce  «  B~'. 

The  second  assumption  Leads  us  to  expect  the  exterior  cold-plasma 
density  to  vary  as  since 

n  ~  JN_  ~  N  ~  for  constant  N. 

Vd.  B^ 

UUe  apply  this  approach  to  the  RMPHED  magnetic  field  operated  in  its 
flat-field  configuration. 
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